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ABSTRACT 

This  report  presents  an  analysis  of  the  impact  errors  of  a  ship-launched  ballistic 
missile  produced  by  an  anomalous  gravity  field.  All  other  error  sources  are  assumed 
to  be  absent.  The  normal  earth  for  which  the  impact  errors  would  be  zero  is 
regarded  as  spherical  and  non-rotating.  The  standard  atmosphere  is  accounted  for. 
Both  navigation  and  in-flight  effects  are  considered,  the  former  under  the  assumption 
that  the  ship  remains  stationary  at  the  launch  point. 

The  treatment  is  largely  deterministic,  with  elementary  statistical  concepts  being 
introduced  to  derive  a  description  of  the  impact  errors  applicable  to  the  case  of  an 
unknown  anomalous  field. 

Some  numerical  statistical  results  are  presented,  in  the  form  of  a  table  for  a 
particular  trajectory  and  as  graphs  for  trajectories  having  various  ranges  and  times  of 
fiight. 


Possible  extensions  of  the  analysis  are  discussed. 
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1.  INTRODUCTION 


This  report  presents  the  analytical  methods  and  results  of  an  investigation  to 
determine  the  extent  of  gravity  data  required  for  a  ship-launched  ballistic  missile  to 
achieve  specified  impact  accuracy.  The  position  and  velocity  of  the  ship  at  missile 
launching  are  assumed  to  be  derived  from  an  inertial  navigation  system  and  the 
trajectory  of  the  missile  is  assumed  to  be  controlled  by  an  inertial  guidance  system. 
Since  the  operation  of  both  of  these  inertial  systems  depends  on  knowledge  of  the 
earth’s  gravity  field,  incomplete  or  erroneous  knowledge  of  this  field  leads  to  both 
navigation  and  in-flight  errors  which  produce  the  resultant  impact  errors. 

To  allow  a  concentration  on  gravity  effects  without  the  complications  that 
would  be  introduced  by  a  simultaneous  consideration  of  other  error  sources,  it  is 
assumed  that  the  ship’s  inertial  navigation  system  and  the  missile’s  inertial  guidance 
system  are  perfect  in  all  respects  except  for  the  omission  of  information  on  the 
anomalous  gravity  field.  More  specifically,  it  is  supposed  that  if  the  field  were 
strictly  normal,  then  this  normal  field  could  be  taken  into  account  completely  and 
no  impact  errors  would  result.  The  impact  errors  considered  are  therefore  those 
which  result  because  the  anomalous  field  is  unknown  or,  if  known,  is  deliberately 
not  accounted  for. 

It  is,  of  course,  impossible  to  make  definitive  statements  about  the  impact 
errors  produced  by  an  unknown  anomalous  field.  The  analysis  therefore  proceeds  as 
though  the  anomalous  field  were  known  but  is  deliberately  not  taken  into  account 
in  navigation  and  guidance.  Formulas  are  then  obtained  which  would  permit  the 
computation  of  the  resulting  impact  error  for  a  specific  set  of  launch  conditions  and 
a  specific  set  of  parameters  describing  the  anomalous  field.  By  averaging  these  results 
over  a  certain  ensemble  of  launch  conditions  a  statistical  description  of  the  impact 
errors  is  obtained.  This  description  is  found  to  have  such  a  form  that  it  does  not 
depend  on  a  detailed  description  of  the  anomalous  field  but  only  on  the  partial 
description  contained  in  the  so-called  “degree  variances”  of  geoidal  heights  or, 
equivalently,  on  the  “autocovariance  function”  of  geoidal  heights.  It  is  this  fact 
which  allows  the  results  to  be  given  a  plausible  interpretation  for  the  case  where 
the  anomalous  field  has  not  been  taken  into  account  because  it  is  not  completely 
known,  since  it  is  not  unreasonable  to  assume  that  the  degree  variances  or  the 
autocovariance  function  may  be  approximately  known  even  though  a  complete 
description  of  the  anomalous  field  is  unavailable. 

The  statistical  results  can  also  be  interpreted  as  describing  the  impact  errors 
that  would  result  if  a  known  part  of  the  anomalous  field  (rather  than  merely  the 
normal  field)  were  taken  into  account  in  navigation  and  guidance  and  the  remainder 


of  the  anomalous  field  were  neglected.  Other  interpretations  of  the  results  will  be 

explained  when  the  results  are  presented. 

In  addition  to  assuming  that  the  navigation  and  guidance  systems  contain 
perfect  components,  it  is  supposed  that  the  ship  remains  stationary  for  a 

considerable  period  (one  or  two-  hours)  prior  to  launching  the  missile.  This  situation 
does  not  seem  unrealistic  and  leads  to  a  simplification  of  the  analysis  in  which 
certain  dynamic  effects  in  the  ship’s  inertial  navigation  system  need  not  be 
considered.  It  may  be  possible  to  avoid  this  restriction  by  an  extension  of  the 

methods  to  be  described. 

It  is  believed  that  the  results  obtained  give  rather  reliable  estimates  of  the 

minimum  amount  of  gravity  data  which  would  be  required  to  achieve  specified 
impact  errors.  It  would  be  expected  that  imperfect  components  in  the  navigation 

and  guidance  systems,  and  a  non-zero  speed  of  the  ship,  would  lead  to  larger 

impact  errors  than  those  derived  here,  but  that  such  a  degradation  could  be 

alleviated  by  the  employment  of  more  extensive  gravity  data. 

Other  analyses  of  this  problem  have  taken  into  account,  the  speed  of  the  ship 
and  the  use  in  the  navigation  system  of  imperfect  components  (particularly  gyros, 
necessitating  a  consideration  of  navigation  resets).  However,  these  treatments  require 
much  more  elaborate  statistical  techniques  whose  validity  is  difficult  to  evaluate.  One 
of  the  objectives  of  the  present  investigation  has  been  to  employ  deterministic 
mathematical  methods  to  the  greatest  extent  possible  and  then  to  introduce 
elementary  statistical  concepts  only  in  the  final  stages. 


II.  THE  ANOMALOUS  FIELD 


The  anomalous  gravity  field  which  is  the  source  of  the  errors  considered  in  this 
report  can  be  described  in  various  ways.  That  which  seems  to  be  most  satisfying 
intuitively,  and  is  chosen  here,  is  the  geoidal  height  function  N(0,X)  of  geodetic 
colatiiude  0  and  longitude  X'  which  gives  the  height  of  the  geoid  above  the 
reference  ellipsoid  of  revolution  corresponding  to  the  normal  gravity  field.  The 
rotation  of  the  earth  introduces  many  complications  into  the  computation  of  impact 
errors  but,  it  is  believed,  does  not  contribute  significantly  to  the  impact  errors 
produced  by  the  anomalous  field.  To  avoid  these  complications  the  rotation  of  the 
earth  is  neglected  in  this  investigation.  This  is  done  by  interpreting  the  given  geoidal 
height  function  N(0,X)  as  defining,  at  geocentric  colatitude  0  and  longitude  X,  the 
height  of  a  non-rotating  geoid  above  a  non-rotating  reference  sphere;  and  then 
replacing  the  normal  gravity  field  corresponding  to  the  rotating  reference  ellipsoid  by 
an  inverse  square  gravitational  field.  With  these  simplifications,  it  is  supposed  that 
the  navigation  and  guidance  systems  account  perfectly  for  the  inverse  square  field 
but  do  not  account  for  the  anomalous  field  corresponding  to  the  geoidal  height 
function. 

Positions  are  defined,  as  illustrated  in  Figure  1,  in  terms  of  spherical 
2  coordinates  r,  0,  X  or  the  corresponding 

a  rectangular  coordinates 


FIGURE  1 


y 


x  =  r  sin  0  cos  X 
y  =  r  sin  0  sin  X  (1 ) 

z  =  r  cos  0  . 


The  origin  is  at  the  center  of  mass  of  the 
earth,  the  z-axis  is  directed  toward  the 
north  pole,  and  the  x-axis  is  in  the 
Greenwich  Meridian. 

The  gravitational  potential  of  the 
earth  will  be  denoted  by  U(r,0,X)  or,  in 


rectangular  coordinates,  by  u(x,y,z).  The 
spherical  coordinate  form  of  the  potential  function  can  be  expanded  in  a  series  of 
solid  spherical  harmonics 
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U(r,0,X) 


a„T(« 


,X)] 


where 


i u  = 

R  = 
a"1  = 

n 


Y™(0,X) 


earth’s  gravitational  constant 

a  reference  length  taken  as  the  radius  of  the  reference  sphere 

complex  coefficients  defining  the  field,  and 

are  the  complex  surface  spherical  harmonics  defined  by 


Y™(0,X)  =  (-1)" 


rsrrr 

(a-  *n)t 

|  4ir 

(n  +  m)! 

P™  (cos  0)eim  x 


(2) 


in  which  PJJ’tz)  is  the  associated  Legendre  function 


Pm(z)  =  - 

n  '  2"n! 


(1  -  z2) 


dn  +  m(z2  -  l)n 


dz 


n  +  m 


(3) 


The  functions  Y™(0,X)  satisfy  the  relations  Y~m(0,X)  =  (- 1  )m  Y™  (0,X)  where  the  bar 

denotes  the  complex  conjugate.  From  this  and  the  fact  that  U  is  real  valued,  it 

follows  that  the  complex  coefficients  must  satisfy  the  conditions  a'm  =  (- 1  )m  ajj1 . 

hi  the  sum  over  n  the  terms  for  n  =  1  have  been  omitted  as  a  consequence  of  the 

origin’s  being  at  the  earth’s  center  of  mass.  The  limits  of  the  summation  indices, 

and  the  indices  themselves,  are  often  omitted  in  what  follows  when  thej^  can  be 

inferred  from  the  context.  It  should  be  admitted  that  the  infinite  sum  £  raises 

n=  2 

questions  of  convergence  in  some  of  the  following  operations.  In  practice,  however, 
the  potential  is  never  known  exactly  and  will  be  approximated  by  a  finite  sum  for 
which  questions  of  convergence  do  not  arise.  The  numerical  values  of  p  and  R 
adopted  here  are 


fj  =  3.986012  X  1014  m3/sec2 
R  =  6.371008608  X  106  m. 
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The  nautical  mile,  used  only  in  describing  the  range  of  a  trajectory,  is  regarded  as 
the  distance  on  a  sphere  of  radius  R  subtended  by  a  central  angle  of  one  minute 
of  arc. 

The  potential  of  the  normal  (inverse  square)  field  is 


v(x,y,z)  =  V(r,0,X)  =  7 


The  disturbing  potential  is 


w(x,y,z)  =  u(x,y,z)  -  v(x,y,z) 


or.  in  spherical  coordinates. 


W(r,0,X)  =  U(r,0,X)  -  V(r,0,X) 


g 

r 


2  <Y-(0,X). 

m  —  n 


(4) 


The  geoid  is  defined  as  the  surface  on  which  the  actual  potential  U(r.0,X)  has 
the  same  value  that  the  normal  potential  V(r,0,X)  has  on  the  surface  of  the 
reference  sphere.  The  equation  of  the  geoid  is  therefore 


u 

R 


r  «  /i 

^  \  n  +n 

1  +  2  (- 

n=2  \  I 

-)  2  a^Y™(0,X) 

■/  m  =  -n 

Let 


r  =  R  +  N  =  R(1  +  ^) 

where  N  =  N(0,X)  is  the  gcoidal  height.  In  terms  of  N  the  equation  of  the  geoid  is 
then 
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Throughout  the  remainder  of  this  report  the  coefficients  ajj’  are  regarded  as  small 
quantities  of  the  first  order  whose  products  and  powers  can  be  neglected.  With  this 
approximation. 


N  =  R  2  2*  a'fY'fffl.X) .  (5) 

If  the  coefficients  a™  in  the  expansion  of  the  potential  are  given,  values  of  N(0,X) 
can  be  computed  from  this  equation.  If,  on  the  other  hand,  values  of  the  geoidal 
height  function  N(0,X)  are  given  over  the  surface  of  the  reference  sphere,  then  the 
coefficients  ojj1  can  be  computed  from 

J(*2 n  p  n 

I  N(e,X)Y?  (0,X)  sin  0  d0  dX  .  (6) 

\=0  Je= o 


The  geoidal  height  function  is  arbitrary  except  for  the  restrictions  that  the  resulting 
values  of  o;®  =  a'1  =  =  a}  =  0.  The  condition  a[j  =  0  implies  that  the  mean 

value  of  N(0,X)  over  the  surface  of  the  sphere  is  zero.  The  three  conditions 
a, 1  =  =  aj  -  0  imply  that  the  centroid  of  the  volume  enclosed  by  the  geoid  is 

at  the  origin  of  the  coordinate  system. 
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III.  EQUATIONS  OF  VARIATION 


A  nominal  missile  is  defined  as  one  which  is  launched  from  a  known  position 
on  the  reference  sphere,  flies  in  the  normal  gravitational  field,  and  has  zero  impact 
error.  The  first  situation  considered  is  that  in  which  a  nominal  missile  is  launched 
from  the  north  pole  at  a  target  in  the  Greenwich  Meridian.  Its  equations  of  motion 
are 


3v 

x°  =  a»<t)+  aii 


-  0  + 


9v 

9y 


o 

9v 

z0  =V,)+ta 


where  av(t),  a„(t),  a,(t)  are  the  components  of  the  non-gravitational  acceleration 
produced  by  thrust  and  aerodynamic  forces,  with  ay(t)  =  0  because  the  nominal 
trajectory  lies  in  the  zx-plane.  The  quantities 


dv 

3v 

3v 

9x 

o’  ^ 

o  ’  9z 

are  the  components  of  the  normal  gravitational  acceleration  at  the  position 
(x0,y0,z0)  of  the  missile.  For  the  trajectories  considered  in  the  present  study,  the 
thrust  is  terminated  above  the  atmosphere  so  that  ax (t)  and  az(t)  are  also  zero  from 
the  end  of  powered  flight  until  reentry  into  the  atmosphere. 

Corresponding  to  this  nominal  missile  we  consider  the  trajectory  of  an  actual 
missile  which  is  launched  from  the  “indicated”  north  pole  (that  is,  from  the 

position  at  which  the  ship’s  inertial  navigation  system,  as  affected  by  the  anomalous 

gravitational  field,  indicates  that  the  ship  is  at  the  north  pole  of  the  reference 

sphere)  and  whose  flight  is  affected  by  the  anomalous  field.  A  simplified  model  of 
the  missile’s  inertial  guidance  system  is  assumed.  In  this  model,  the  missile  is 
assumed  to  be  controlled  during  powered  flight  in  such  a  manner  that  its 

components  of  non-gravitational  acceleration  are  the  same  functions  of  time  as  are 
the  corresponding  components  for  the  nominal  missile.  The  equations  of  motion  for 
the  actual  missile  during  powered  flight  are  therefore 
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r 


3u 

x  =  a  (t)  +  — 
x  ox 


3u 

y  =  °+  — 

3y 


3u 

z  =  a  (t)  +  — 
z  3z 


From  the  termination  of  powered  flight  until  reentry  these  equations  still  apply  with 
ax(t)  =  az(t)  =  0,  as  in  the  case  of  the  nominal  missile.  From  reentry  until  impact, 
a  (t)  and  az(t)  for  the  actual  missile  are  not  the  same  functions  as  for  the  nominal 
missile,  and  ay(t)  is  not  equal  to  zero. 

Considering  the  portion  of  the  flight  from  launch  to  reentry,  let  £  =  x-xQ, 
V  -  y  _  y0 ,  f  =  z  -  zQ  be  the  components  of  displacement  of  the  actual  missile 
from  the  nominal.  The  differential  equations  satisfied  by  £,  tj,  are  then 


l  = 

3u 

3v 

_ 

3v 

3v 

+ 

3w 

3x 

3x 

0 

3x 

3x 

0 

3x 

77  = 

3u 

3v 

3v 

3v 

+ 

3w 

’i 

3y 

3y 

0 

3y 

3y 

0 

3y 

r  = 

3u 

3v 

3v 

3v 

+ 

3w 

3z 

3z 

0 

3z 

3z 

0 

3z 

If  v  and  w  are  expanded  as  power  series  in  £,  rj,  J  and  terms  involving  powers 
and  products  of  these  quantities  are  neglected,  the  differential  equations  become 
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In  the  case  of  the  normal  gravitational  field,  the  coefficients  aJJ1  are  zero  and  the 
quantities  £,  77,  J  would  be  zero.  This  implies  that  the  products  etc.,  which 

would  arise  from  the  presence  of  the  second  partial  derivatives  of  w  in  the  right 
members  are  small  quantities  of  order  higher  than  the  first  and  may  therefore  be 
neglected.  With  this  simplification  the  equations  become 


32  v 
*  ~  3^2 

32  v 

t  +  - 

0  dyd* 

32  v 

0  77  3z3x 

3w 

f  +  — 

0  9x 

32v 

3x3y 

,  a2v 

S  +  7T 
0  ay 

32v 
v  + 

0  3zdy 

0  8> 

••  32v 

3x3z 

32v 

£  + - 

0  dydz 

32v 

T?+  ry 

0  3z  1 

3w 

T 

)  3z 

As  was  stated  earlier,  the  non-gravitational  acceleration  of  the  actual  missile 
during  reentry  is  not  the  same  as  that  of  the  nominal  missile  and  so,  for  this 
portion  of  the  flight,  the  differential  equations  for  £,  77,  f  must  be  altered.  This 
change  requires  a  more  detailed  consideration  of  the  trajectory.  For  the  terminal 
portion  of  the  flight  of  the  actual  missile  it  is  assumed  that  the  aerodynamic  force 
consists  only  of  a  drag  force  D  which  acts  in  a  direction  opposite  to  the  velocity. 
The  equations  of  motion  are  then  found  to  be 

D  x  3u 

x  = - +  — 

m  v  3x 

D  y  3u 

y  - -  +  — 

m  v  3y 

D  z  3u 

z  - +  ~  • 

m  v  3z 


(Note  that  in  these  equations  v  represents  the  speed  of  the  missile,  not  the  normal 
gravitational  potential.  In  some  of  the  following  equations  v  occurs  with  both 
meanings,  but  the  proper  interpretation  will  be  clear  from  the  context.)  The 
corresponding  equations  for  the  nominal  missile  are 
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_  Dq  Xo  +  9v 
m  v0  3x  0 


Dp  yp.  +  ^ 


_  Do  zo  +  dv 
m  vQ  dz  o 


If  £,  77,  f  are  defined  as  before,  they  satisfy  the  differential  equations 


••  _  1  /  Dx  D0x0\  3u  3v 

m  \  v  v„  )  dx  dx 


_  J_  /Dy  D0y0\  +  3u  3v 
m  \v  v0  /  3y  3y 

■■  _  J_  /Dz  DqZq\  +  du  _  dv 
m  \v  vQ  /  3z  dz 


The  differences  (3u/3x)  -  (3v/3x|0 ) ,  etc.,  are  expanded  as  before  to  obtain  the 
equations 


1  Dx  Dnx, 


32v  32v 

3x2  0  ^  3y3: 


3y3x  0  3z3x  0  dx 


dw 

0  0 


m  \  v 


1  / Dy  Ppyp  \  +  d2v  ^  {  92v 


1  /  DZ  Da  Zr 


dydz  0 


32v  3w 

3z3y  0  3y 


32v  3w 

3z2  /+  dz  o 


Now  (Dx/v)  -  (D0x0/v0)  =  5(Dx/v) ,  etc.,  where  the  expression  on  the  right  denotes 
the  increment  of  Dx/v  in  passing  from  a  point  on  the  nominal  trajectory  to  the 
corresponding  (same  time)  point  on  the  actual  trajectory.  Hence 


1  /Dx\  32  v  32v 
m  \v/  3x2  Q  3y3x 


.  3w 

?  +  r- 


1  c  /Dy\  '  32v  32v  32v  3w 

m  \v  /  3x3y  Q  3y2  0  3z3y  Q  3y 


1  ''  z\  32v  32v  32v 

—  5  —  I  +  -  £  +  -  17  +  — m* 

m  \  v/  3x3z  Q  3y3z  0  dz2, 


32v  .  3w 

TT  ^  +  "a — 

3z  v  0  3z  0 


Furthermore 


6(v7  =  ?6*  +  *6(t) 

D0  .  .  /  Dy 

=  ~  £ +  x0  M 


where  small  quantities  of  higher  order  have  been  neglected.  With  similar 
approximations  in  the  other  two  equations,  the  differential  equations  become 


3x2  Q  3y9x  0  3z3x  0  mv0  m  \v/  3x 


32  v  +  3^_v  32v  Dq 

3x3y  0  3y2  J1  3z3y  0  mv0 


v0  m  \  v/  3y  0 
1  .•  Zn  .  /D\  3w  I 


3x3z  0  3y3z  0  3z2  0 


n+r-T  — f-  “  +T- 


The  drag  force  D  is  assumed  to  be  of  the  form 


D  =  -pv2SCD 


where  p  is  the  atmospheric  density,  CD  is  the  drag  coefficient  of  the  missile,  and  S 
is  its  cross  section  area.  Hence,  D/v  =  (S/2 )pvCD  and 


6 


(vCD  6p  +  pCD  5v  +  pv6CD)  . 


Now  CD  =  Cd(M)  is  a  function  of  the  Mach  number  M  =  v/a  in  which  a  is  the 
speed  of  sound.  Hence 


6Cd  =  Cp(M)6M  =  ^(fiv-  M5a) 


where  =  Cp(M)  is  the  derivative  of  CD  with  respect  to  its  argument  M.  Then 
pv 5Cd  =  pMCp(5v-  M6a)  and 

5 


Also  p  =  p(h)  and  a  =  a(h)  are  functions  of  the  height  h  so  that  6 p  =  p'(h)5h  and 
5a  =  a'(h)5h,  where  the  primes  here  denote  derivatives  with  respect  to  h.  With  these 
substitutions 


=  ^[vCD5P+^CD+MCD)5v"PM2CD5a]  • 


s 

-  [(vCDp'-  pM2Cpa')6h  +  p(CD  +  MCjj  )6y]  , 


which  can  be  abbreviated  as  5 


mP  p,  ,  mQ  c  .  .  .  . 

—  6h  +  - ov  in  which 

v  v 


p  =  ^(vCDp'-pM2C^a> 
Q  "  ^  (CD  +  MCp  )pv  . 


(8) 


The  height  h  =  r  -  R  where  r  is  the  distance  of  the  missile  from  the  center  of  the 
earth  and  R  is  the  radius  of  the  reference  sphere.  Hence  6h  =  6r  and,  since 
r2.  =  x2  +y2  +z2, 


6h‘  =  -  6x  +  -  5y  +  -  5z  . 
r  r  r 


Similarly,  since  v2  =  x2  +y2  +z2, 


5v  =  -  5x  +  -  5y  +  -5z  . 

V  V  V 


Hence 


i'fe)  ■ 


P/x  „  y  .  z  \  Q  /x  „ .  y  r .  z  . 
-  I  -  8x  +  -5y+-6zl  +  — [  -  5x+-6y+  -  5z 
v  r  r  r  v  v  v  v 


P  Q  .  ■  .. 

=  —  (x£  +  yr?  +  zf)  +  — 2  (x£  +  yr?  +  zf ) . 
rv  v4 


Since  one  of  the  small  quantities  appears  in  each  term  of  this  expression, 

the  quantities  r,x,y,z,v,x,y,z,P,Q,  can  be  evaluated  on  the  nominal  trajectory  for 
which  y  =  y  =  0,  and  therefore 


J_ 

m 


-  (x£  +  z$:)+  ■?  (x^  +  zf). 
rv  v^ 


The  differences  between  the  position  and  velocity  components  of  the  actua1  and 
nominal  missiles  are  now  accounted  for  by  the  quantities  p,f;  and  the 

subscript  zero  can  be  omitted  everywhere  with  the  understanding  that  all  quantities 
are  to  be  evaluated  on  the  nominal  trajectory.  The  equations  for  then  become 
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Now  let  x/r  =  sin  0  and  z/r  =  cos0,  where  0  is  the  angle  from  the  z-axis  to 
the  radius  vector  of  the  nominal  missile;  and  let  x/v  =  sin  <p  and  z/v  =  cos^>,  where 
yf>  is  the  angle  between  the  z-axis  and  the  velocity  vector  of  the  nominal  missile. 
Considering  the  normal  potential  v  =  #i/>/x2  +  y2  +  z2  ,  the  partial  derivatives  of  v 
appearing  in  the  differential  equations  are  found  to  be 


92v  h  , 

^2  =  -^(l-3sin'20) 


92v  ji_ 
dy2  ~  ~  r* 


92v  _  n 

~  _pr  (1  -  3  cos2  0) 


92v  =  92v 

9x3y  9y9x 

92v  _  92v 

9y9z  9z9y 


92v  _  92v 

9z9x  9x9z 


~r  sin  0  cos  0  . 
r 
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The  differential  equations  for  £,  tj,  f  can  therefore  be  written  in  the  form 


'  ■  ‘  [5 


(1-3  sin2  0)  +  P  sin 


.  ,  .  1 t  ,  f3M  . 

in  0  sin  y  £  +  — r  si 

J  Lr 


—r  sin  0  cos  0  -  P  cos  0  sin  ^  f 


in  JJ| 


■[™+QsinJ*’) 


+  Q  sin2  \fi  ||j  -  Q  sin  >p  cos  f  +  — - 

ox 


/i  D  .  9w 

T]  ~  3  V - T?+  — 

r  mv  oy 


T 

f  =  +  —  sin  0  cos  0  -  P  sin  0  cos  if 
-  Q  sin 


M?0- 


3  cos2  0 )  +  P  cos  0  COS  i p  K 


in  cos  ip  £  *  —  +  Q  cos2  t  +  — 
l_mv  v  *_|5  az 


It  is  convenient  to  write  these  differential  equations  in  matrix  form  by  putting 
£  =  x,,  tj  =  x2,  f  =  x3,  £  =  x4,  tj  =  x5,  f  =  x,.  This  permits  them  to  be  written 
as  the  single  differential  equation 


dx 

dt 


F(t)x  +  G(t) 


(9) 


where  x  and  G(t)  are  the  columns 


0 

X1 

0 

0 

X2 

3w 

X3 

G(t)  = 

dx 

X4 

3w 

X5 

3y 

9w 

_X6_ 

dz 

and  F(t)  is  the  6  X  6  matrix 


Fit)  = 


0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

Pr  ( 1  -  3  sin2  0 ) 

Mr  sin  6  cos  0 

Tn  1 

1 

P  sin  0  sin  ^ 

0 

U 

-  P  cos  8  sin 

-  [^*0^ 

0 

0 

iL 

0 

0 

_D 

0  -  Q  sin  $  cos  y 


pj  sin  8  cos  0  -  ( 1  -  3  cos2  8) 

—  P  sin  0  cos  +  P  cos  8  cos  ysj 


-  Q  sin  cos  a  0  -  I  —  +  Q  cos2u? 

|_mv  ’j 


The  matrix  F(t)  is  determined  by  the  nominal  trajectory  and  the  column 
G(t)  by  the  partial  derivatives  of  the  disturbing  potential  w  evaluated  on  the 
nominal  trajectory.  A  comparison  of  this  matrix  differential  equation  with  the 
differential  equations  (7)  for  the  portion  of  the  flight  prior  to  reentry  shows  that 
the  matrix  differential  equation  may  be  employed  for  the  entire  flight  if  P,  Q,  and 
D  are  taken  as  zero  for  the  portion  of  the  flight  prior  to  reentry. 


The  solution  of  the  differential  equation  (9)  is 


x  ( t )  =  X(t 


,t„)x(t0)  +  /  x(t, 
J  tn 


r)G(r)dr , 


where  x(tQ)  is  the  value  of  x  at  the  time  tQ  of  launch  and  where  the  6X6  matrix 
X(t.r)  is  the  solution  of  the  homogeneous  matrix  differential  equation 


subject  to  the  initial  condition 


dX 

dt 


=  F(t)X 


X(t,t)  =1  =  6X6  identity  matrix. 


The  time  tQ  of  launch  will  in  practice  be  taken  as  zero  but  is  retained  as  tQ  in 
some  of  the  equations  for  purposes  of  clarity.  We  will  be  interested  in  the  value  of 
x(t)  at  the  time  t  =  T  of  impact  of  the  nominal  missile,  that  is,  in 


x(T)  =  X(T,t0)x(t0)  + 


X(T,t)G(t)dt . 


GO) 


It  can  be  shown  that  X(T,t)  =  YT(t,T),  where  the  superscript  T  denotes  matrix 
transpose  and  Y(t,T)  is  the  solution  of  the  “adjoint”  matrix  differentia!  equation 


dY 

dt 


-FT(t)Y 


(ID 


subject  to  the  “initial”  condition 


Y(T,T)  =  I  . 


For  use  in  Equation  (10),  X(T,t)  is  required  not  only  for  t  =  tQ  but,  to  evaluate 
the  integral,  for  all  values  of  t  in  the  range  tQ  <  t  <  T.  To  obtain  these  matrices, 
the  homogeneous  differential  equation  which  is  actually  solved  is  that  for  Y(t,T), 
starting  at  t  =  T  with  the  identity  matrix  and  integrating  backwards  in  time  to 
t  =  tQ.  This  gives  Y(t,T)  over  the  range  T  >  t  >  t0  and  hence  X(T,t)  =  YT(t,T) 
over  the  range  tQ  <  t  <  T.  An  examination  of  the  form  of  F(t)  shows  that  X(T,t) 
must  have  zero  elements  in  the  positions  indicated  in  the  expression 


X(T,t)  = 


X,, 

0 

X1  3 

X1  4 

0 

X1  6 

0 

^2 

0 

0 

X2  5 

0 

X31 

0 

X3  3 

X34 

0 

X36 

X41 

0 

X4  3 

X44 

0 

X46 

0 

X5  2 

0 

0 

X5  5 

0 

X61 

0 

X6  3 

X64 

0 

X66 

♦ 


♦ 


m 
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IV.  IMPACT  ERRORS 


The  column  x(T)  is  formed  from  the  differences  £(T),  i?(T),  f(T),  £(T),  r?(T), 
f(T)  between  the  position  and  velocity  components  of  the  actual  missile  and  those 
of  the  nominal  missile  at  the  instant  of  impact  of  the  nominal  missile.  Figure  2 
illustrates  the  nominal  trajectory 


FIGURE  2 


in  the  zx-plane.  The  nominal  missile  strikes  the  target  on  the  reference  sphere  at 
range  angle  0(T).  At  this  instant  the  actual  missile  is  displaced  from  the  nominal 
missile  by  a  component  £(T)  in  the  x-direction,  a  component  f(T)  in  the  z-direction, 
and  a  component  i?(T)  (not  shown  in  the  diagram)  in  the  y-dircction.  During  the 
additional  time  to  impact  of  the  actual  missile  it  will  move  from  this  position. 
Since  only  first  order  effects  are  being  considered  in  the  analysis,  the  velocity 
of  the  actual  missile  during  this  small  time  interval  is  regarded  as  being  equal  to  the 
impact  velocity  of  the  nominal  missile.  This  velocity  v.  forms  an  angle  \p  with  the 
downward  vertical  at  the  target. 

Since  _v  lies  in  the  zx-plane,  the  actual  missile  will  have  no  motion  in  the 
y -direction  during  the  small  time  increment  under  consideration  and  the  cross-range 
component  of  impact  error  will  be  5C  =  r?(T),  taken  as  positive  to  the  left  when 
viewed  from  the  launch  point.  At  the  nominal  impact  time,  the  actual  missile  is 
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above  the  target  a  distance  £(T)  sin  0(T)  +  f(T)  cos  0(T)  and  is  down  range  from  the 
target  a  distance  £(T)  cos  0(T)  -  f(T)  sin  0(T).  Since  the  vertical  component  of 
velocity  of  the  missile  is  v  cos  \j/,  it  will  impact  at  a  time 
[£(T)  sin  0(T)  +  f(T)  cos  0(T)]/(v  cos  1/1)  later  than  the  nominal  missile.  During  this 
time  its  down-range  component  of  velocity  is  v  sin  \p  and  hence  it  will  travel  a 
distance  |£(T)  sin  0(T)  +  f(T)  cos  0(T)1  tan  \p  in  the  down-range  direction.  Adding  this 
to  the  down-range  distance  from  the  target  at  the  nominal  impact  time,  the 
down-range  component  of  impact  error  is 


6D  =  £(T)  cos  0(T)  -  f(T)  sin  0(T) 

+  [£(T)  sin  0(T)  +  f(T)  cos  0(T)]  tan 

=  a«T)  +  bf(T) 


where 


a  =  tan  \p  sin  0(T)  +  cos  0(T) 
b  =  tan  \p  cos0(T)-  sin  0(T)  . 

The  angles  0(T)  and  i p  can  be  expressed  in  terms  of  the  x  and  z  components  of 
position  and  velocity  of  the  nominal  missile  at  time  t  =  T  and  a  and  b  expressed 
in  the  forms 


a 

b 


R 

+  — : - r  z 

xx  +  zz 

R  . 

- : - r  x. 

xx  +  zz 


(12) 


Because  the  anomalous  field  and  the  geoidal  height  have  been  defined  in  terms  of 
the  complex  surface  spherical  harmonics  Y^(0,X)  with  the  complex  coefficients  ckJJ1  , 
it  becomes  convenient  to  combine  the  impact  error  components  into  a  complex 
impact  error 


p  -  5D  +  i8C 


=  a£(T)  +  irj(T)  +  bf(T) 


=  [a  i  b] 


tnT 

t?(T) 

f(T) 


From  this  equation  it  is  seen  that  only  the  First  three  elements  of  x(T)  are 
required  and  hence  only  the  first  three  rows  of  the  matrix  X(T,t).  Furthermore, 
under  the  assumption  that  the  ship  has  remained  stationary  for  a  considerable  period 
before  the  missile  is  launched,  and  the  further  assumption  that  the  ship’s  inertial 
navigation  system  operates  in  a  “damped”  mode  involving  the  use  of  non-inertial 
horizontal  velocity  data  and  height  (or  depth)  data,  the  velocity  outputs  of  the 
navigation  system  will  be  zero  so  that  £(t0)  =  r?( tQ)  =  f(t0)  =  0.  Therefore,  in  the 
product  X(T,t0)  x(tQ)  in  Equation  (10),  only  the  first  three  columns  of  X(T,t)  are 
needed.  In  the  integrand  function  X(T,t)G(t)  of  the  same  equation,  the  first  three 
elements  of  G(t)  are  zero  so  only  the  last  three  columns  of  X(T,t)  are  required. 
Making  use  of  these  facts  it  is  found  that  p  can  be  written  in  the  form 


p  =  jaX, ,  +  bXj,  X22  aXl3+bX„] 


J  (aXM+bX34  X2S  aX,6+bX36J 


% 

-f°_ 

dw 

dx 

dw 

'  dy 
dw 


dt 


(13) 


where,  now,  £0  =  £(t0),  etc.  The  matrix  elements  X^  in  the  first  row  stand  for 
Xjj(T,tfl ),  which  are  constants  for  a  particular  nominal  trajectory.  The  matrix 
elements  which  appear  in  the  row  under  the  integral  sign  stand  for  Xjj (T,t),  which 
are  functions  of  time  along  the  nominal  trajectory. 


The  two  rows  (1X3  matrices)  appearing  in  Equation  (13)  involve  a,b,  and  the 
matrix  elements  X^  and  depend  only  on  the  nominal  trajectory.  The  two  columns 
depend  on  the  anomalous  field.  The  first  column;  whose  elements  are  £0,  i7?0,  ?0; 
represents  the  initial  position  errors  due  to  the  anomalous  field  and  is  independent 
of  the  trajectory.  When  multiplied  by  the  trajectory-dependent  row 
[aX,  j  +  bX3 ,  X22aX,  3  +  bX33  ] ,  the  contribution  to  p  of  the  initial  position  errors 
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(i.e.,  ship's  navigation  errors)  is  obtained.  The  second  column;  whose  elements  are 
3w/3w,  idw/3y,  3w/3z;  depends  on  the  anomalous  field  (since  w  is  the 
disturbing  potential)  and  also  depends  on  the  nominal  trajectory,  since  the  partial 
derivatives  are  to  be  evaluated  on  this  trajectory.  This  column  is  then  multiplied  by 
the  trajectory-dependent  row  [aX14+bXJ4  X2J  aX16+bX36)  and  the  product  is 
integrated.  This  integral,  being  independent  of  the  initial  position  errors,  represents 
the  in-flight  contribution  to  the  complex  impact  error  p.  In  the  next  two  sections 
the  navigation  and  in-flight  contributions  to  p  are  evaluated  separately. 


V.  NAVIGATION  CONTRIBUTION  TO  IMPACT  ERROR 


The  launch  point  of  the  nominal  missile  is  the  north  pole  of  the  reference 
sphere  and  the  ship  attempts  to  bring  itself  to  this  position  by  using  the  outputs  of 
its  inertial  navigation  system.  The  initial  position  errors  £0,  are  the 

component  displacements  of  the  ship  from  the  north  pole  of  the  reference  sphere 
when  the  erroneous  outputs  of  the  navigation  system  indicate  that  it  is  at  the  north 
pole  of  the  reference  sphere.  Since  the  navigation  system  contains  no  information  on 
the  geoidal  height  function,  the  actual  position  of  the  ship  will  be  on  the  geoid 
instead  of  on  the  reference  sphere.  Since  the  navigation  system  also  contains  no 
information  on  the  vertical  deflections  and  operates  in  a  damped  mode,  the  vertical 
indication  of  the  navigation  system  will  be  the  plumb  line  vertical  which  is  normal 
to  the  geoid.  It  has  also  been  assumed  that  the  inertial  navigator  employs  perfect 
gyroscopes.  It  can  therefore  be  supposed  that  it  possesses  axes  in  the  x,y,z 
directions  whose  orientation  has  been  accurately  established  prior  to  the  mission  by 
an  in-port  alignment  procedure  and  that  this  orientation  is  maintained  throughout 
the  mission  by  the  perfect  gyros.  Since  the  ship  is  believed  to  be  at  the  north  pole 
of  the  reference  sphere,  its  actual  position  will  therefore  be  at  such  a  place  on  the 
geoid  that  the  normal  to  the  geoid  is  parallel  to  the  z-axis. 

From  Equations  (1)  the  parametric  equations  of  the  geoid  are 


x  =  (R  +  N)  sin  0  cos  X 

y  =  (R  +  N)  sin  0  sin  X  (14) 

z  =  (R  +  N)  cos  0 


where  N  =  N(0,X)  is  given  in  terms  of  0,  X,  and  the  o^1  by  Equation  (5).  If  0  and 
X  are  given  small  increments  d0  and  dX,  the  increment  of  z  will  be 


dz 


■cos  0  -  (R  +  N)  sin  6  d0  + 


in  0  c 


9N 

ax 


cos0dX. 


If  the  increments  d0  and  dX  are  from  such  values  of  0  and  X  that  the  normal  to 
the  geoid  is  parallel  to  the  z-axis,  then  dz  =  0  and  hence  the  coefficients  of  d0 
and  dX  arc  zero,  or 
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—  cos  6  =  (R  +  N)  sin  6 
bd 


3N  . 

—  cos  6  =  0  . 

ax 


(15) 


These  are  two  simultaneous  equations  for  6  and  X.  If  they  can  be  solved  and  the 
solution  values  substituted  in  Equations  (14),  the  coordinates  of  the  ship  will  be 
obtained  and  the  initial  position  errors  will  be  (since  the  intended  coordinates  are 
0,0. R) 


il0  =  (R  +  N)  sin  0  cos  X 

r)0  =  (R  +  N)  sin  d  sin  X  (16) 

?0  =  (R  +  N)  cos  0  -  R  . 


It  is  possible  that  Equations  (15)  have  many  solutions  for  d  and  X  and  that 
some  of  these  solutions  might  have  large  values  of  6  and  represent  ship’s  positions 
far  from  the  north  pole.  This  statement  corresponds  to  the  fact  that  it  is  possible 
conceptually,  although  it  is  probably  not  actually  the  case,  that  there  may  be 
distinct  points  on  the  geoid  having  the  same  astronomic  coordinates,  i.e.,  distinct 
points  where  the  plumb  lines  are  parallel.  However,  if  the  ajj1  are  all  zero  the  geoid 
coincides  with  the  reference  sphere  and  the  only  solution  (except  for  the  south 
pole)  is  0  =  0.  Hence,  if  the  aJJ1  are  small  there  must  be  a  solution  having  a  small 
value  of  0  and  this  solution  will  be  sought  as  the  one  appropriate  to  the  problem. 

Since  0  vanishes  with  the  aJJ1  and  powers  and  products  of  these  small 
quantities  are  being  neglected,  the  approximations  sin  6  =  6  and  cos  6  -  1  can  be 
made  in  Equations  (15)  to  obtain 


9N 

b6 

8N 

ax 


=  (R  +  N)0 
=  0. 


(17) 
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The  functions  N,  3N/30,  3N/3X  appearing  in  these  equations  can  be  expanded 
as  power  series  in  0  (X  being  regarded  as  a  parameter)  of  the  forms 


N  =  Nlo+^  • 


3N 

3NI 

32N 

30 

h 

_ o 

1  ^ 

1  fO 

3»2l 

3N 

3N 

82  NL 

303X|O 

3X 

3X0  + 

in  which  terms  involving  powers  of  0  higher  than  the  first  have  been  neglected.  The 
bar  on  the  right  indicates  an  evaluation  at  0  =  0  which  may,  however,  depend  on 
X.  Substituting  these  expansions  in  Equations  (17)  and  again  omitting  terms  involving 
02 ,  one  obtains 


3N  32N 

an  +  9  (R  +  N|o)0 


30 10  30  2  |0 


3N|  d2N  „ 

—  + - 0-0 

3X|0  303Xjo 


To  obtain  N  and  its  partial  derivatives  evaluated  at  0  =  0,  Equation  (5)  is  first 
rewritten  with  the  substitution  of  the  expressions  (2)  for  the  Y^”  (0 ,X): 


N  =  R  2  I«”(-ir  m);P!?(cos0)eimX 

n  m  n  N  4?r  (n  +  m)!  n 


Differentiation  with  respect  to  0  yields 


3N  m  m  |2n+I  (n-m)!  dP™(cos0).. 

—  =  R2  Za"(-ir  \| - - - !L~ ^ - e  x 

30  n  m  n  '  4tt  (n  +  m)!  d0 
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and,  using 


dPlT  (cos  9 )  1  A1  1 

sin0 - - —  =  -P"1  '(cos®)-  ~(n  +  m)(n  -  m  +  l)P™_1(cos0) ,  (20) 

d(cos0)  2  n  2  n 


this  becomes 


BN  R  ^  m  ,vm  |2n  +  l  (n  -  m)!  r„m . .  n 

90  "  '  2n  m“n("  U  J  4tt  (n  +  m)!  P*"  ~  (n  +  m)(n  -  m  +  1)P™ -|]  eim  * 


In  this  and  some  later  equations  the  argument  of  the  P™  is  always  to  be 
understood  as  cos0  when  it  is  not  exhibited  explicitly.  Differentiating  (19)  with 
respect  to  A  gives 


=  iR2  2  ajj1  (- 1  )m  m 


|2n  1 

N  At, 


111  lniinl!pmeimx 

4?r  (n+  m)!  n 


and  the  second  partial  derivative 


92N  =  R, 
909A  ”  1 2‘i 


=  _i2?  Z7ff[C+1-(n  +  mKn-m  +  1)C_1]eir 


To  evaluate  these  expressions  at  0  =  0,  use  is  made  of  the  relation,  from  (2), 

(-  l)m  Y?  (0,A)e"im  x 
P"1  (cos  0)  =  °  ;  ==- 

n  pn  +  1  (n  -  m)! 

N  47T  (n  +  m)! 


and  the  fact  that 


V?(0,X)=  ^6” 


where  8™  is  the  Kronecker  delta.  With  the  abbreviation  A,,  = 
results  are 


n(n  +  l)(2n  +  1) 
4jt 


the 


■Nl0 


R2 

n  * 


2n+  1 


4?r 


3N 

30' 


Anan 

n  n 


i  _ 


i\ 


2  A  Ofn 

n 


f 


3N 


=  o 


_a l 

303X1, 


J]  .  R 

[e,x  2  Aa1  +  e 

ix  2  Ancr  *1 

o  =‘'2 

L  n  n  n 

n  11  nJ 

(21) 


The  second  partial  derivative  32N/302|O  has  not  been  evaluated  because  it  is 
multiplied  by  0  in  the  first  of  Equations  (18)  and  this  product  contains  no  first 
order  terms  in  the  ajj1 .  This  product  and  also  the  product  N|Q0  in  the  same 
equation  are  therefore  omitted.  With  these  approximations,  the  substitution  of  the 
values  (21)  in  Equations  (18)  yields 


R  r  ~i 

-  —  elA  2  A  a2  -  e_,x2  A  a" 1 
2  L  n  n  n  n  n  n  J 

R  r 

-i  —  elX  2  Ana2  +  e~,x2  AnoC 1 
2  I  "  n  n  n  J 


=  R0 

0  =  0. 


(22) 


We  are  seeking  a  value  of  0  =£  0  which  satisfies  these  equations  and  hence 


elXZ  AnaJ  +  e_,x2  A_a  1  =0 

n  n  n  n  n  n 


(23) 


from  which  X  could  be  determined.  With  this  value  of  X  the  first  of  Equations  (22) 
would  give  the  value  of  0.  These  two  values  could  be  substituted  in  the  first  two 
of  Equations  (16)  to  obtain  and  r?0.  However,  these  equations  can  be  combined 
into  the  form 


£0  +  ir?0  =  (R  +  N)  sin  0(cos  X  +  i  sin  X)  =  (R  +  N)  sin  0  elX 
or,  to  the  desired  degree  of  approximation, 

£0  +  if?0  =  R  0  eix  . 

Using  the  first  of  Equations  (22), 

f0 A»“!  - e‘k 

But  X  satisfies  (23),  which  can  be  written 

e2ix  2  A  od  =  -  l  A norl 

n  n  n  n  n  n 

and  hence 

$0  +i?70  =  RSA,,^1 
=  -R2  AaJ  . 

_  n  n 
n 

The  conjugate  of  this  expression  is 


*0"  "»0  =  “R£  AX  ' 


The  last  two  equations  can  be  solved  for  £0  and  ir/0  to  obtain 


*0  =  -  ^  AnK  +  =  -  I?  AM  "  “n1  > 


ij?0  =  +  ff  An(“n  -  =  +  f  ?  AnK  +  “n  ‘  > 


2  n 


(24) 


We  also  need  the  component  f0  of  initial  position  error  which,  from  (16),  is 


=  (R  +  N)  cos  6  -  R 


or,  approximately, 


•  ni„  ■  *5^jF«e 


(25) 


For  use  in  later  developments  it  is  desirable  to  express  Equations  (24)  and  (25)  as 
double  sums  (over  m  as  well  as  n)  by  introducing  the  Kronecker  deltas.  Recalling 
the  definition  of  A^,  these  equations  then  become 


In  = 


-s-.) 

+  +s->) 

_  |2n  +  1 


5™ 
4jt  0 


(26) 


We  now  return  to  Equation  (13)  and  define  trajectory-dependent  constants 
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(27) 


gj  -  3Xj  j  (T,tg )  +  bXj  j  (Tjtg ) 

~  X22(T,to) 

g3  =  aX j  3  (T,t0 )  +  bXj  3  (T,t0 ) . 

Using  these  constants,  the  navigation  contribution  to  the  complex  impact  error  can 
be  written  first  as 


p(Nav)  =  gj£0  +  g2if?0  +  g3?0 


and  then,  by  introducing  the  values  (26),  as 


p(Nav) 


av)  =  ?  J 

n  m  3 


2n  +  1 
4jt 


>/n(n  +  1 ) 


-  Rgr^-r-^  -  5-) 


(28) 


+  Rg2  (6y  +  s™, )  +  Rg3  sg1  . 
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VI.  IN-FLIGHT  CONTRIBUTION  TO  IMPACT  ERROR 


This  contribution  to  the  complex  impact  error  arises  from  the  integral 
expression  in  Equation  (13).  It  has  been  explained  how  the  1  X  3  matrix  in  the 
integrand  can  be  obtained,  as  a  function  of  t,  by  integration  of  the  adjoint 
Equation  (11).  We  now  consider  the  column  involving  the  partial  derivatives  of  the 
disturbing  potential  w.  These  partial  derivatives  are  related  to  the  partial  derivatives 
of  the  spherical  coordinate  form  W  of  the  disturbing  potential  by  the  transformation 


3w 

7>e_  3x  af 

~3W 

9x 

9x  9x  9x 

90 

dw 

90  9X  9r 

9W 

dy 

9y  3y  9y 

9X 

dw 

90  9X  9r 

9W 

3z_ 

_3z  3z  9z_ 

_9r_ 

The  3X3  matrix  is  evaluated  by  making  use  of  the  relations  (1)  between 
rectangular  and  spherical  coordinates  and  remembering  that,  because  the  partial 
derivatives  are  to  be  evaluated  on  the  nominal  trajectory  for  which  X  -  0,  the 
partial  derivatives  in  the  3X3  matrix  can  be  evaluated  for  X  =  0.  Also  introducing 
the  factor  i  we  obtain 


9w 

_ 

1  9W 

— 

cos  0  0  +sin  0 

—  - 

3x 

r  30 

.  9w 

i  9W 

1  - 

= 

0  1  0 

-  ■  — 

3y 

r  sin  0  9X 

9w 

9W 

-  sin  0  0  cos  0 

— 

3z 

3r 

_  — 

—  __ 

The  integrand  involved  in  the  complex  impact  error  therefore  becomes 


j 

i 

< 

i 

i 

cos  0 

0 

+sin  0 

r 

9W 

90 

[aX|4+bX34  X25  aX16+bX36l 

0 

1 
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9W 
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i 
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0 
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By  introducing  the  functions 


/j  (t)  =  (aX,  4  +  bX34 )  cos  0  -  (aX,  6  +  bX36 )  sin  0 

/,(!)=  X„  (29) 

/3(t)  =  (aX,4  +  bX34)sin0  +  (aX]6  +  bX36)cos0 


the  integrand  may  be  written  as 


1  3W  i 

/,(»)-  —  +/2(t)— - 
1  r  30  L  r  sin  0 


3W 

3X 


+  /3(  0 


aw 

3r 


The  functions  /£ ( t )  depend  on  t  not  only  through  the  factors  sin  0  and  cos0  but 
also  because  the  matrix  elements  X^  =  Xjj(T,t). 

The  partial  derivatives  of  W  are  found  by  using  the  expression  (4)  for  W(r,0,X) 
and  the  definitions  of  Y™(0,X)  and  PJJ’fcos 0)  in  (2)  and  (3).  In  addition, 
formula  (20)  and  the  recursion  formula 


P™ + 1  (cos  0) -  2m  cot  0P™  (cos  0)  +  (n  +  m)(n  -  m  +  OP™-1  (cos  0)  =  0 


are  required.  After  forming  the  partial  derivatives,  X  is  set  equal  to  zero,  since  the 
derivatives  are  to  be  evaluated  on  the  nominal  trajectory.  The  results  of  these 
operations  are 
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r  30 
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n 


Z  (-  1  )m 


2n  +  1 
4tr 


(n  -  m)! 
(n  +  m)! 


(m  cot  0  •  P">  -  P™  + 1 

n  n 
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i  3W 
r  sin  0  3X 


m  n  1  4ff 


(n  -  m)! 
(n  +  m)! 


m  esc  0  •  P™ 


3W 

3r 


-  2  (n  +  1 ) 

r  n 
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2  o%  (-  1  )m 
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2n  +  1 
4  ?r 


<"  ~  m>!  pm 

(n  +  m)!  n 
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Substituting  these  expansions  in  the  integrand  function  and  interchanging  the  orders 
of  integration  and  summation,  the  in-flight  contribution  to  the  complex  impact  error 
becomes 


P(I-  F)  =  2 


With  the  substitutions 


this  becomes 

p( I  -  F) 


2<(-l)mpRn 


2n  +  1  (n  -  m)! 
47r  (n  +  m)! 


[/,  (m  cot  dP™  -  P^  + 


/2  m  esc  0  P™  -  /3  (n  +  1  )P™J  dt . 


J/, (t)(m  cot  0  •  -  P™  +  1 )-  /3(t)(n  +  1)P™J  dt 


p  .  f  ’  — 

rn  +  ‘ 


(30) 


J?  =  I  /2(t)  m  esc  6  ■  P™dt 

*0 


=  2  2a"'(-I)m 

n  m  n 


_  2n  +  1  (n  -  m)!  _ 

“R  >|— 


(31) 


VII.  THE  TOTAL  COMPLEX  IMPACT  ERROR 


The  total  impact  error  p  in  Equation  (13)  is  the  sum  of  the  navigation 
contribution  p(Nav)  in  (28)  and  the  in-flight  contribution  p(I  -  F)  in  (31).  This  sum 
can  be  written  in  the  form 


P-SK  c 

n  m 


m 

n 


(32) 


where 


2n  + 

cm  =  \z - 

n  =  '  4n 


>Jn(  n  +  1 )  Vn(n  +  1 ) 

-  Rg,  ~  (5™  -  )  +  Rg2  ~  — (5?  +  5m, )  +  Rg3  5™ 


2 

+  (- 1  r 


pRnJ  —  "  Jn)l  • 

N  (n  +  m)!Vn  n/J 


(33) 


For  use  in  some  later  operations  it  should  be  noted  that  the  c™  are  real. 

To  summarize  the  development  thus  far,  suppose  that  the  nominal  missile  is 
launched  from  the  north  pole  of  the  reference  sphere  and  strikes  a  target  in  the 
Greenwich  Meridian.  The  radius  vector  r  and  range  angle  6  of  this  missile  are 
known  func  ions  r  =  r(t)  and  0  =  0(t)  of  the  time,  with  impact  occurring  at  time 
t  =  T.  By  integrating  the  adjoint  Equation  (11)  corresponding  to  this  trajectory,  the 
constants  g( ,  g2.  g3  defined  by  Equations  (27)  and  the  functions  /,  (t),  /2(t),  /3(t) 
defined  by  Equations  (29)  are  determined.  The  integrals  IJJ1  and  JJJ1  defined  by 
Equations  (30)  cun  therefore  be  evaluated  and  the  quantities  cj]1  computed.  The 
complex  impact  error  of  the  actual  missile  is  then  given  by  Equation  (32),  where 
the  aJJ1  are  the  coefficients  in  the  expansion  (5)  of  the  geoidal  height  function.  The 
impact  error  thus  derived  includes  the  effect  of  the  initial  position  error  resulting 
from  the  action  of  the  anomalous  field  on  the  ship’s  inertial  navigation  system  and 
the  in-flight  effect  due  to  the  action  of  the  anomalous  field  on  the  flight  of  the 
missile. 


VIII.  ARBITRARY  LAUNCH  POSITION  AND  TARGET  AZIMUTH 


It  will  next  be  explained  how  the  previous  results  can  be  applied  to  the 
determination  of  the  impact  error  of  a  missile  which  is  launched  from  an  arbitrary 
position  on  the  earth  toward  a  target  at  an  arbitrary  azimuth.  The  principle  which 
is  involved  in  tiiis  extension  of  the  results  depends  on  the  fact  that  in  the  case 
already  considered  the  nominal  launch  position  was  required  to  be  at  the  north  pole 
of  the  reference  sphere  and  the  target  in  the  meridian  of  Greenwich  only  because 
the  expansion  of  the  geoidal  height  function  was  made  in  the  terms  of  the 
functions  Yni  (0,X)  of  colaiitude  0  (measured  from  the  north  pole)  and  longitude  X 
(measured  from  the  Greenwich  Meridian),  and  the  coefficients  cdj1  involved  in  the 
expansion  were  based  on  the  use  of  these  particular  functions.  If  we  wish  to 
consider  a  launch  point  other  than  the  north  pole,  we  can  regard  the  launch  point 
as  the  fictitious  “north  pole”  of  a  new  coordinate  system  and  the  plane  of  the 
nominal  trajectory  as  the  fictitious  “Greenwich  Meridian”  of  the  new  coordinate 
system.  In  this  new  coordinate  system,  the  geoidal  height  would  naturally  be 
expressed  as  a  function  N'  =  N'(O'.X')  of  “colatitude”  O'  and  “longitude”  X',  and 
expanded  in  terms  of  the  functions  Y™(0',X').  The  coefficients  in  this  expansion 
would  be  certain  constants  instead  of  the  c^1  in  the  original  expansion.  If  these 
constants  can  be  determined,  the  complex  impact  error  for  the  new  situation  can  be 
computed,  by  analogy  with  Equation  (32).  from  the  equation 


P 


=  2  Z  0"c” 

n  m  "  n 


(34) 


where  the  constants  cjj1  have  the  same  values  that  they  have  for  the  nominal 
trajectory  originating  at  the  north  pole  and  lying  in  the  Greenwich  Meridian. 

To  carry  out  this  idea,  let  a  represent  the  longitude  of  the  new  launch 
point,  (3  its  colatitude,  and  7  the  azimuth  of  the  target,  measured  CCW  from  the 
south.  The  relation  of  the  new  coordinate  system  (x\  y',  z'  or  O',  X')  to  that 
originally  used  (x.y.z  or  0,X)  is  indicated  in  Figure  3,  while  Figure  4  shows  how 
the  “colatitude'  O'  and  “longitude”  X'  of  an  arbitrary  point  on  the  reference  sphere 
are  related  to  the  true  colatitude  0  and  true  longitude  X  of  the  same  point.  The 
surface  spherical  harmonics  YJJ,(0,,X')  in  the  new  coordinate  system  are  related  to 
those  in  the  old  system  (M.  E.  Rose,  Elementary  Theory  of  Angular  Momentum, 
Equation  (4.28a))  by  the  equation 
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RELATION  BETWEEN  X.Y.Z  AXES  ANDX'.Y'.Z'  AXES 


z 
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RELATION  BETWEEN  $,  X  COORDINATES  AND 
z  9[  X1  COORDINATES 


FIGURE  4 


Y"m(°')=  k£n  dL^)yJ;(0,x), 


(35) 


where  the  functions  D£m  have  a  simple  dependence  on  a  and  7  but  a  rather 
complicated  dependence  on  0.  The  actual  form  of  these  functions  will  not,  however, 
be  needed.  The  inverse  relationship  is 

YfJ1  (0,X)  =  2  DjJm  (-7,-/3, -a)Y^(0',X') 

which,  by  virtue  of  Rose’s  Equation  (4.21),  may  be  written 

Y™  (0.X)  =  2  D^k(aJJ.7)  Yj \(0'X)  . 

The  geoidal  height  expansion  (5)  can  therefore  by  written  as 
N(0.X)  =  R2  JDMaF)Yk(0') 

n  m  k  m  n 

- R  f  Y:<9'.vj]. 

This  is  an  expansion  of  the  geoidal  height  in  the  form 

N'(O')  -  R  2  2  0™  Y™  (O'X) 

n  m 

and  hence 


C  =  0™(«M)  =  2^D|Jm(a,(3(7)  • 

The  complex  impact  error  for  the  case  now  under  consideration  can  therefore  be 
written 


=  2  2  2  c™ 

n  m  k 


<%DnkmMy). 


(36) 
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IX.  STATISTICAL  CONSIDERATIONS 


Up  to  this  point,  the  development  of  an  expression  for  the  impact  error  has 
been  entirely  deterministic.  We  might  be  given  at  the  outset  a  particular  nominal 
trajectory  and  could  compute  the  values  of  the  c“  given  by  Equation  (33)  without 
regard  to  the  launch  position  and  target  azimuth.  We  would  also  be  given  the  values 
of  the  ajj1  in  the  expansion  of  the  geoidal  height.  A  particular  launch  point  would 
then  be  specified  by  the  values  of  a  and  /3  and  a  particular  target  azimuth  by  the 
value  of  7.  Using  the  well-defined  functions  D£m(cii,  0,  7)  whose  explicit  form  has 
not,  however,  been  displayed,  the  coefficients  (3™  could  be  computed.  Finally,  the 
complex  impact  error  would  be  computed  from  Equation  (34). 

Although  all  of  these  computations  might  actually  be  carried  out,  we  are  in 
practice  often  interested  in  obtaining  only  a  statistical  description  of  the  impact 
errors  for  some  ensemble  of  launch  conditions.  The  statistical  description  which 
would  be  most  useful  would  be  to  consider  a  particular  target  and  to  describe 
statistically  the  impact  errors  of  missiles  launched  toward  this  target  from  various 
points  of  some  limited  launch  area.  If  the  earth’s  gravitational  field  were  known, 
but  had  deliberately  been  neglected,  this  problem  could  be  solved  by  a  Monte  Carlo 
method  in  which  a  numerical  statistical  analysis  would  be  made  of  the  impact 
errors,  computed  as  just  outlined,  of  missiles  launched  from  a  large  number  of 
randomly  selected  points  within  the  launch  area.  Such  a  computation  would  be 
extremely  time-consuming  as  a  result  of  the  large  number  of  nominal  trajectories 
(many  sets  of  the  cJJ1 )  and  the  large  number  of  launch  points  (many  sets  of  the 
)  that  would  have  to  be  considered  to  obtain  an  adequate  sample  of  the  impact 
errors.  Furthermore,  the  results  obtained  would  apply  to  the  case  where  the  known 
field  had  deliberately  been  neglected,  a  situation  in  which  the  impact  errors  could 
have  been  avoided  by  taking  that  field  into  account.  By  considering  a  different 
ensemble  of  launch  conditions,  however,  it  becomes  possible  to  obtain  a  statistical 
description  of  the  impact  errors  by  a  method  which  is  largely  analytical  and  involves 
numerical  computation  to  a  much  lesser  extent.  Of  greater  importance,  this  method 
is  applicable  to  the  situation  where  the  anomalous  field  has  not  been  taken  into 
account  because  it  is  not  completely  known. 

The  situation  considered  is  one  in  which  a  particular  nominal  trajectory  is 
chosen;  defined,  for  example,  by  range  and  time  of  flight.  Then  it  is  supposed  that 
missiles  are  launched  from  positions  uniformly  distributed  over  the  surface  of  the 
earth  and  that,  for  each  such  launch  point,  missiles  are  launched  at  targets  having  a 
range  from  the  launch  point  equal  to  the  range  of  the  nominal  trajectory  and  being 
uniformly  distributed  in  azimuth.  A  statistical  description  of  the  impact  errors  in 
such  a  situation  is  the  goal  of  the  remainder  of  this  report. 
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The  first  statistic  to  be  computed  is  the  expected  or  mean  value  p  of  the 
complex  impact  error  p(otJ},y).  This  has  the  value 


P  = 


8ir2 


/"  /'  / 

Ja  =  0  3  =  0  Jy-Q 


2  it 


p(otJi,y)  sin  /3  dy  d/3  da . 


We  will  frequently  be  concerned  with  integrals  of  this  form  and  adopt  the 
abbreviation 


/2it  r  n  *  2it 

I  I  f(a#,y)  sin  0  d7  d/3  da  =  J/(a,0,7)df2 
0  •'0  =  0  Jy  =  0 


When  /(a,/3,7)  =  1  the  integral  becomes  jdSl  =  8tt2. 
Hence 


P  =  ^2  Jp(aJ3»7)dn 


From  Equation  (36)  this  becomes 


P  = 


■  s?  5  s 


To  evaluate  the  integral,  use  is  made  of  Rose’s  Equation  (4.60): 


Sir2 


(37) 
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Now,  from  Equation  (35), 


Yg(0\X')  =  S  D®0(a^,7)Y *(0,X)  =  D° 0 «*,/), 7 )Y°(0 ,X) . 

But  Yj|(0,X)  =  \l  1  /4jt  Pg(cos0)  =  sj  1/4*  and  also  Y[j(0',X')  =  v/ I/4n  ,  so  that 
Dq0(O!,(3,7)  =  1.  Hence,  putting  n'  =  m'  =  k'  =  0  in  (37), 


jD£m(a,/3,7)dft 


0—2 

^77  W5  -  sir^Dsys; 


and  the  expression  for  p  becomes 


p  =  -^  £  L  2  cj"a!'8ff26?5J!1 =  0 

8tt2  n  -  2  m  k  n  "  0  0  0 

because^  n  >  2.  Since  p  -  0,  its  real  and  imaginary  parts  are  also  zero  and 
6D  =  6C  =  0.  Hence,  the  mean  values  of  the  down-range  and  cross-range 
components  of  the  impact  error  are  both  zero. 


The  remaining  statistics  to  be  computed  are  the  variance  8D2  of  the 
down-range  impact  errors,  the  variance  SC2  of  the  cross-range  impact  errors,  and  the 
covariance  5D  •  5C.  In  terms  of  p  and  its  conjugate  p,  the  down-range  and 
cross-range  components  of  impact  error  are 


5D  -  +  ~(p  +  p) 
5C  =  -  ~2  (p  -  p) . 


Hence 


5D2  = 


1 

+  “(p  +  p)(p  +  p) 

4 
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5C2  =  -7 (p  -  p)(p  -  p) 

4 

5D  •  8C  =  -  -<p  +  p)(p  -  p) 

4 

so  that 

6D5  =  +  7— yf(p  +  p)(p +  p)dft 

5C2  =  ~  TTT  ib> "  p)(p  "  P)dn 

32?r  J 

5D^5C  =  -  ^yj (p  +  p)(p  -  p)dft  . 


The  reduction  of  these  expressions  to  the  desired  forms  follows  the  general  pattern 
used  in  reducing  the  expression  for  p.  The  computations  are  tedious  but  offer  no 
special  difficulties.  The  properties  of  the  D£m(a,0,7)  that  must  be  made  use  of  are 
Equation  (37)  and 


jDknni(a,/?,7)  Dnk:m,(a^,7)dS2  =  (-l)k'"m' 


8tt2 
2n+  1 


n  -  m 


This  last  expression  is  obtained  by  replacing  m1  and  k'  in  Equation  (37)  with  -m' 
and  -k'  and  then  applying  Rose’s  Equation  (4.22)  to  substitute 


Dn'k- _m.(a£,y)  =  (-  l)k  m  Dk'm’,Of’0’'*) 


-k',-m 


As  a  result  of  these  operations  it  is  found  that 


6D*  =  2  ak 
n  n  n 

6C5  =  I  bk 
n  n  n 


5D  *  5C  =  0  , 


(38) 
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where  the  real  positive  quantities 


4  it 


+  n  _ 

2  cPa”* 

nc _  n  ^  ^ 


(39) 


are  the  so-called  “degree  variances”  of  the  geoidal  height.  The  quantities  an  and  bn 
depend  only  on  the  nominal  trajectory  and  are  given  in  terms  of  the  cm  of 
Equation  (33)  by 


a 


n 


2tt  1 
R2  2n  +  1 

2ff  1 

R2  2n  + 1 


[c” 

[c”  -  (-I)mc;”]  . 


Using  the  expression  for  the  c™  derived  earlier,  these  sums  can  be  expressed  as 
sums  over  non-negative  values  of  m  with  the  results 


an  *  +  7  n(n  +  1  )gj  +2f*R"~l  (g3 1°  +  gj  I„ ) 


+  H 


2r2(ii-1)  |(]0)2  +  2  £  (P)21 

L  n  m  =  i  (n  +  m)!  n  J 


(40) 


bn  =  -n(n+  l)g2  +2^Rn_1(g2J‘) 


+  ^RJ<n-„  2  |  1SJJS21  y:)! 

m  =  i  (n  +  m)!  n 


To  compute  the  an  and  bn  for  a  specific  trajectory  it  is  necessary  first  to 
obtain  r  and  6  as  functions  of  the  time;  then  to  integrate  the  adjoint  Equation  (11) 
to  obtain  the  constants  gp  g3  and  the  functions  /j  ft),  /2(t),  /3(t);  and.  finally, 
to  evaluate  numerically  the  integrals  and  JJJ1  in  (30).  The  most  difficult  part  of 
this  task  is  to  compute  the  associated  Legendre  functions  PJJ1  (cos  <9 )  for  the  values 
of  6  =  0(t)  occurring  along  the  trajectory.  Some  preliminary  computations  of  these 
functions  for  large  values  of  n  (and  m,  since  0  <  m  <  n)  suggested  that  the 
computations  could  be  performed  more  accurately  in  terms  of  the  normalized 
functions  PJ^cosfl)  (the  bar  here  denotes  normalization,  not  complex  conjugation, 
the  latter  process  not  being  involved  in  any  of  the  subsequent  operations)  defined 
by 
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pm  =  r \n±±  -  m)!  nm 
n  **  2  (n+m)!  n 


When  these  functions  are  introduced  it  is  found  that  Equations  (40)  can  be  written 
in  the  forms 


an  =  g\  +  -n(n+  l)g2  +  2MRn~1(g3I°  +gjNn(n+  1)  I*  ) 


+  m2r2(i-1)  j|,0)2  +2  £  fl®)^ 

bn  -  ~n(n  +  l)g2  +  2/LtRn' *(g2  \/n(n  +  DJ*) 
+  /x2R2(n-l)  £  £  ^  (J®)2j, 


where  the  “normalized”  integrals  I™  and  J™  are  defined  by 
I™  =N|  (n  -  m)!/(n  +  m)!  IJJ1  and  JJJ1  =  yj  (n  -  m)!/(n  +  m)T  J™  and  may  be  written 
in  the  forms 


=>2^Tl  /  ^2  [/,  («(»  «*  e  '  Pn  -  >/(n-  ™  +  1  >  pr l) 

^‘o  -  f3(t)(n  +  l)PjJdt 

■C  =>J^Ti  /  PTT5  /2<t>  m  CSC 0  -  ^dt. 

‘o 

In  computing  the  integrand  functions  the  factors  cot  0  and  esc  0  would  be 
expected  to  cause  difficulty  with  the  small  values  of  6  occurring  near  the  launch 
point,  except_  for  m  =  0  when  the  terms  involving  these  factors  are  absent.  But  for 
m^Ojhe  P™(cos0)  carry  sin  6  as  a  factor  and  the  products  cot  6  •  PJJ*  (cos  0)  and 
esc  0  •  P™(cos0)  remain  finite  as  0  approaches  zero.  To  make  use  of  this  fact, 
another  set  of  functions 


Q™(cos0)  =  esc  0  •  P™(cos0) 
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is  defined  (these  functions  should  not  be  confused  with  the  Legendre  functions  of 
the  second  kind).  The  integrals  then  become 

C  [/j(t)(mcos0  •<£  ->17n-  m)(n  +  m+l)P^T1) 

**°  -/3(t)(n+l)P^]dt 

The  values  of  the  P"1  are  needed  for  0  <  m  <  n  and  the  values  of  the  Qm 

n 

are  needed  for_l  <  m  <  n.  For  each  value  of  0  occurring  along  the  trajectory,  the 
values  of  the  PjJ1  and  the  Q™  were  computed  by  the  following  recursion  schemes: 


For  values  of  0  <  n  <  N  and  0  <  m  <  n, 

*5  -J-  *  m  I-6 

n  =  2,3.-,N  i  JRTTJtsrny*..^- 

nrn  —— 

n  =  l^.'^.N  P"  =  — —  sin0*PJJ"} 


n 


2,3,-",N 


pn-  1 
n 


>/2nTT 

n 


cos  6  •  P""}  + 


n  -  1  2n  +  1 
n  >|2n  -  2 


sin  0 


m  =  1  :  n  =  3,4. -'’.N 

m  =  2  :  n  =  4,5,*  ".N 


•  •  • 

m  =  N-  2:n  =  N 


2n  -  1  (2n+  l)(n-  m) 

n  >  (2n-  l)(n  +  m) 


cos  6  •  P™_ , 


n-  1  (2n+  l)(n -  m)(n-  m-  1)  — — 

n  >  (2n  -  3)(n  +  m)(n  +  m  -  1)  n~2 


2n  -  1 


n 


m 


J 


2n+  1 


(2n-  l)(n  +  m)(n  +  m-  1) 


sin  6 


pm  - 1 
rn-  1 
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For  values  of  1  <  n  <  N  and  1  <  m  <  n, 


n  =  2,3,*  *\N 


n  =  2,3,***,N 

n  =  3,4,-  •  \N 

m  =  2  :  n  =  4,5,“,,N 

in  =  3  :  n  =  5,6,',#,N 


ff-Jf 


Qj  =  -^ycos0 


qT=  *z±  (2n  +  "  jj  cos  g  »  qT 

n  n-  1  J(2n-  l)(n  +  1)  n_1 


"  — 
n-  1  J 


(2n+  l)(n  -  l)(n  -  2) 
(2n  -  3)(n  +  1  )n 


Q1  . 

vn-  i 


Q"  '  J^r sin  9 '  ^ 


\n-  1 


J  2n  +  1 


cos  0  .  Qn:  1  +  1_L  sin  0  .  Qn^T 


— —  2n  - 

cr  = 


-  1  R 2nT 

n  J(2n- 

n-  1  f 
n  W 


l)(n  -  m) 


l)(n  +  m) 


cos  6  •  QJ^_  j 


(2n  +  1  )(n  -  m)(n  -  m  -  1 ) 

(2n  -  3Xn  +  m)(n  +  m  -  1)  v"“2 


Q?.; 


m  =  N  -  2  :  n  =  N 

I 


2n-  1 


m 


J 


2n+  1 


(2n  -  1  )(n  +  m)(n  +  m  -  1) 


sin  6  ■  Q™_ 


m  -  1 
1 


I 


* 

! 
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X. 


TRAJECTORY  COMPUTATIONS 


For  the  trajectory  computations  it  was  assumed  that  the  missile  consisted  of 


powered 

stages,  designated  as  Stages  1  and  2,  and  a  reentry 

body  .designated  as 

3.  The 

gross  characteristics 

of  the  three  stages 

are  listed  in 

the  following  table 

Stage 

Initial  Mass 

Rate  of  Mass 

Thrust 

Cross-Section 

Decrease 

Area 

(kg) 

(kg/sec) 

(newtons) 

<m2) 

1 

73754.0 

795. 

2035200. 

6.132 

2 

9330.0 

105. 

286230. 

2.775 

3 

72.1 

0. 

0. 

0.1 17 

The  first  stage  was  assumed  to  burn  for  a  fixed  time  of  75  sec  and  the  second 
stage  for  a  variable  additional  time  T2  from  0  to  65  sec.  To  avoid  considering 
details  of  the  guidance  system,  it  was  assumed  that  the  trajectories  are  zero  lift  in 
which  the  only  aerodynamic  force  is  the  drag  acting  in  a  direction  opposite  to  the 
velocity  and  the  thrust  (for  the  first  two  stages)  acting  in  the  sane  direction  as  the 
velocity.  The  vertical  velocity  at  launch  (tQ  =  0)  was  taken  as  20  m/sec  and  the 
horizontal  velocity  at  launch,  VH ,  was  given  a  value  in  the  range  0-4  m/st\  By 
varying  the  initial  horizontal  velocity  and  the  burning  time  of  the  second  stage,  a 
family  of  trajectories  having  various  ranges  and  total  times  of  (light  is  obtained. 

Drag  data  for  the  three  stages  were  given  initially  in  the  form  of  tables  of  urag 
coefficient  CD  vs.  Mach  number  M.  These  values  were  fitted  with  Walton-Shanks 
formulas  of  the  form 


CD  =  (1  +  J)A(X)  +  (1  -  J)B(X) 
A(X)  =  Aq+AjX  +  AjX2 
B(X)  =  B„  +  B,X 


X 

f(l  -0)X2  +  0 


M2  -  a 
M2  +  a 


The  constants  have  the  following  values. 


Stage 


2 


3 


+8.8212-01 
+  1.0618-02 
+  1.7087-01 
-3.2897-02 
-1.4217-02 
+2.7938-02 
-2.0284-03 


+  1.0808-00 
+2.5982-02 
+2.0935-01 
-9.6134-02 
+2.1966-02 
+  1.0241-02 
- 1 .2650-02 


+8.1741-01 
+2.2353-05 
+2.3078-01 
-3.0271-02 
-1.3220-01 
+1.5429-01 
+  1.7857-02 


The  derivative  required  in  integrating  the  adjoint  equation,  is  obtained  by 

differentiating  the  Walton-Shanks  formula. 

The  atmospheric  density  and  speed  of  sound  as  functions  of  geopotential 
altitude  were  taken  from  the  tables  of  the  U.  S.  Standard  Atmosphere  with  the 
following  modifications.  For  geometric  altitudes  greater  than  90,000  meters,  the 
Standard  Atmosphere  is  defined  in  a  different  and  more  complicated  manner  than  it 
is  for  lower  altitudes.  To  allow  a  uniform  mode  of  computation  to  be  used  at  all 
altitudes,  this  portion  of  the  Standard  Atmosphere  was  altered  slightly.  However,  the 
density  at  these  altitudes  is  so  small  that  the  change  has  a  negligible  effect  on  the 
results.  Another  change  made  in  the  Standard  Atmosphere  was  to  replace  the 
polygonal  function  defining  the  molecular  scale  temperature  TM  as  a  function  of 
geopotential  altitude  by  an  eighth  degree  polynomial  giving  1/TM  as  a  function  of 
geopotential  altitude.  With  these  changes,  the  derivatives  p'(h)  and  a'(h)  required  in 
integrating  the  adjoint  equation  could  be  obtained  analytically.  The  atmospheric 
density  was  taken  as  zero  at  geometric  altitudes  >  150000  meters. 

The  trajectory  computations  were  made  with  a  fourth  order  Runge-Kutta 
integration  scheme.  The  time  increments  were  taken  as  follows:  0.25  sec  from 
t  =  0  to  t  =  5  sec,  1 .25  sec  from  5  sec  to  75  sec,  5  sec  from  75  sec  to  an 
altitude  of  approximately  150,000  meters,  and  0.25  sec  from  that  point  until 
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impact.  One  exception  to  this  scheme  occurs  at  the  termination  of  second  stage 
thrust,  where  two  smaller  intervals  whose  sum  is  5  sec  were  used  to  allow  thrust 
termination  to  occur  at  an  arbitrarily  chosen  time. 

Using  a  set  of  arbitrarily  chosen  “even”  values  of  T2  and  VH ,  a  family  of 
trajectories  was  computed  having  the  ranges  and  times  of  flight  plotted  in  Figure  5. 
Using  numerical  values  of  partial  derivatives  obtained  from  these  trajectories,  values 
of  T2  and  VH  were  derived  by  an  iteration  process  to  produce  the  final  45 
trajectories  listed  in  the  following  table.  In  this  table  T  denotes  the  total 
time-of-flight  in  seconds  and  R  the  range  in  nautical  miles. 


T  R  T2  V„ 


0 

28.68013141 

0.00000000 

500 

29.75426850 

0.55712446 

1000 

32.78836295 

1.1  1086563 

1500 

37.29593791 

1.65948124 

2000 

42.67621520 

2.20437465 

2500 

48.40687876 

2.75286094 

3000 

54.23841118 

3.32571427 

0 

43.16535414 

0.00000000 

500 

43.58334483 

0.41185108 

1000 

44.79389359 

0.82380098 

1500 

46.67771033 

1.23624173 

2000 

49.06795304 

1.65010040 

2500 

51.78537344 

2.06705416 

3000 

54.67051088 

2.48996986 

3500 

57.61244696 

2.92502477 

4000 

60.59935459 

3.3  8334173 

0 

52.21427363 

0.00000000 

500 

52.40753281 

0.33187979 

1000 

52.97353510 

0.66433233 

1500 

53.87333100 

0.99799642 

2000 

55.04924451 

1.33362982 

2500 

56.43370819 

1.67215207 

3000 

57.958351 15 

2.01467201 

3500 

59.56213362 

2.36260283 

4000 

61.19814676 

2.71812983 

4500 

62.84295251 

3.08608459 

5000 

64.52198318 

3.47026982 
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RANGE  AND  TIME  OF  FLIGHT  VS  T2  AND  VH 
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FIGURE  5 


T 

2240 


2560 


R 

T2 

VH 

0 

58.22928795 

0.00000000 

500 

58.33024376 

0.28248605 

1000 

58.62760481 

0.56551324 

1500 

59.10555530 

0.84963652 

2000 

59.73994425 

1.13543461 

2500 

60.50114277 

1.42351394 

3000 

61.35721867 

1.71450737 

3500 

62.27702102 

2.00907542 

4000 

63.23288008 

2.30793586 

4500 

64.20281645 

2.61201070 

0 

62.44253642 

0.00000000 

500 

62.50011875 

0.24938435 

1000 

62.67025943 

0.499208!  1 

1500 

62.94540583 

0.74991005 

2000 

63.3138121! 

1.00192618 

2500 

63.76065880 

1.25568501 

3000 

64.26937279 

1.51159996 

3500 

64.82298548 

1.77006024 

For  all  of  these  trajectories,  the  thrust  is  terminated  at  an  altitude  greater  than 
150,000  meters,  justifying  the  assumption  that  the  non-gravitational  acceleration  is 
zero  from  the  time  of  thrust  termination  to  reentry.  For  each  trajectory  the 
position  and  velocity  components  at  the  time  of  impact  were  used  in  formulas  (12) 
to  obtain  the  values  of  a  and  b  which  were  required  later  to  compute  the  constants 
gj ,  g2,  g3  of  Equations  (27)  and  the  functions  (t),  f2( t),  f3  (t)  defined  by 
Equations  (29). 

To  allow  the  integration  of  the  adjoint  equation  it  would,  in  principle,  be 
necessary  to  retain  from  the  integration  of  the  nominal  trajectory  only  the  values  of 
r(t)  and  0(t)  for  tQ  <  t  <  T,  but  it  is  more  efficient  to  retain  the  values  of  r, 
sin0.  cos0  (0  itself  is  not  actually  required),  sin  «p,  cos<p,  D,  m.  and  v.  In  addition, 
although  not  actually  required  during  the  integration  of  the  nominal  trajectory,  the 
values  of  P  and  Q  are  computed  from  Equations  (8)  and  retained  for  use  during 
the  integration  of  the  adjoint  equation. 

The  adjoint  equation  is  solved  by  a  Runge-Kutta  integration  (backwards  in 
time)  using  an  integration  step  of  0.5  sec  from  t  =  T  to  a  time  near  that  of 
reentry  and  a  step  size  of  10  sec  from  there  to  the  time  of  launch.  This  integration 
gives  Y(t,T)  for  T  >  t  >  tQ  and  hence  X(T,t)  =  YT(t,T)  for  tQ  <  t  <  T.  As  a 
check  on  the  solution  of  the  adjoint  equation,  it  was  solved  for  the  trajectory 
T  =  1920  sec  and  R  =  3000  nautical  miles  and  the  final  matrix  X(T,t0)  printed 
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out.  Then  the  nominal  trajectory  program  was  rerun  and  the  values  of  the 

non-gravitational  acceleration  components  stored  in  the  permanent  file.  Trajectories 
were  then  run  with  varied  initial  conditions  where,  instead  of  computing  the 

non-gravitational  acceleration  components,  they  were  read  from  the  permanent  file 
for  the  times  prior  to  reentry  and  computed  anew  from  reentry  to  the  time  of 
nominal  impact.  From  the  differences  between  the  position  and  velocity  components 
at  the  nominal  time  of  impact  for  the  varied  trajectories  and  the  nominal  trajectory, 
it  was  possible  to  compute  numerical  values  of  16  of  the  20  non-zero  elements  of 
X(T,  t0).  The  remaining  four  non-zero  elements  could  not  be  obtained  by  this 

method  because  the  trajectory  program  did  not  allow  out-of-plane  variations  of  the 
initial  conditions.  The  values  of  the  16  non-zero  elements  of  X(T,t0)  agreed  with 
those  obtained  from  the  solution  of  the  adjoint  equation  to  four  or  more  significant 
figures. 

When  this  method  of  checking  was  first  attempted,  the  polygonal  TM  function 
of  the  Standard  Atmosphere  was  in  use  and  the  two  methods  gave  results  agreeing 

to  only  one  to  three  significant  figures.  It  was  for  this  reason  that  the  polynomial 

approximation  to  1/TM  was  adopted,  since  this  eliminated  the  discontinuities  in  the 
derivatives  occurring  with  the  polygonal  function.  A  similar  difficulty  might  have 

been  experienced  by  employing  the  tables  of  drag  coefficient  originally  available,  but 
the  use  of  the  Walton-Shanks  representation  of  the  drag  coefficient  at  the  outset 
avoided  the  problem. 

in  the  integration  of  the  adjoint  equation,  36  elements  of  the  X(t,tQ)  matrix 
are  determined  (16  of  which  are  zero),  but  not  all  of  these  are  required  to 

compute  the  an  and  bn.  The  elements  actually  required  are  those  appearing  in  the 

functions  /,(t),  /2(t),  /3(t)  defined  by  Equations  (29)  and  those  (for  t  =  tQ  only) 
appearing  in  the  constants  of  Equations  (27).  These  values  were  stored  in  the 

permanent  file  for  later  use  in  the  computation  of  the  an  and  bn.  The  values  of  r, 

sin  0.  cos  0  which  had  been  developed  during  the  integration  of  the  nominal 
trajectory  and  used  during  the  integration  of  the  adjoint  equation  were  also  put  in 
the  permanent  file  along  with  /j(t),  /2( t),  /3(t). 

The  major  paid  of  the_computation  of  the  an  and  bn  consists  of  the  evaluation 
of  the  integrals  I™  and  J™ .  This  was  done  by  Simpson’s  Rule,  using  a  10-sec 

interval.  In  this  stage  of  the  computation  it  was  convenient  to  replace  r,  the  radius 
vector  of  the  missile,  by  the  normalized  value  r/R,  with  R  the  radius  of  the 
reference  sphere;  to  take  R  =  1  in  the  formulas  (41)  for  an  and  bn;  and  to  express 
p  in  these  equations  in  units  of  (reference  sphere  radius)3 /sec2.  The  values  of 

/,  (t),/2(t),/3(t),  r,  sin  0,  cos0  required  to  evaluate  the  integrals  were  read  from  the 
permanent  file  created  during  the  integration  of  the  adjoint  equation.  A  subroutine 
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of  the  program  for  computing  the  an  and  bn  employed  the  recursion  formulas 
which  have  been  presented  earlier  for  the  computation  of  the  P™  and  the  . 

To  allow  a  break-down  of  the  total  errors  into  navigation  and  in-flight  errors, 
an  and  bn  in  Equations  (41)  are  written  in  the  forms 


an  = 


aln  +a2n  +a3n  +a4n 


where 


*1  n 


A 


“2  n 


3n 


4  n 


-  n(n  +  1  )g2 


=  2pRn-i(g3lJ+g1^rr)iI) 


=  /i2R2(n'*) 


R)2+2  2(C)2 

V  n/  m  =  }V  n  / 


bn  =  b2n  +b3n  +b4n 


where 

b2„  =  -n<n  +  l)gJ 

b,„  =  jsrrnij) 

b4n  =  (c)]. 

An  examination  of  the  origins  of  these  terms  shows  that  if  the  effect  of  in-flight 
errors  only  is  to  be  computed,  then  an  should  be  taken  equal  to  a4n  and  bn  equal 
to  b4n.  If  the  effect  of  navigation  errors  only  is  to  be  computed,  then  an  should 
be  taken  equal  to  aln+a2n  and  bn  equal  to  b2n.  The  effect  of  navigation  errors 
only  can  be  further  broken  down  into  the  effect  of  horizontal  errors  only 
(an  =  a2n  and  bn  =  b2n)  and  the  effect  of  height  errors  only  (an  =  aln  and 
bn  =  0).  It  will  be  noted  that  aln,  a2n,  a4n,  b2n,  b4n  are  necessarily  positive 
whereas  aJn  and  b3n  are  not.  JThese  last  two  values  represent  the  effects  on  the 
total  mean  square  impact  errors  5D2  and  6C^  of  correlations  between  the  navigation 
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and  in-flight  effects.  If  a3n  and  b3n  are  zero,  the  total  effect  is  the  statistical  sum 
of  the  separate  navigation  and  in-flight  effects.  If  a3n  or  b3n  is  positive,  the 
navigation  and  in-flight  effects  (for  the  corresponding  component  of  impact  error) 
tend  to  reinforce  one  another  and  if  negative  tend  to  cancel  one  another. 

To  allow  all  of  these  effects  to  be  evaluated  separately,  if  desired,  the  program 
for  computing  the  an  and  bn  for  each  trajectory  stores  on  magnetic  tape  the  value 
of  aln  (actually  independent  of  n)  and  for  each  value  of  n,  2  <  n  <  174,  the 
values  of  a-,  a,„,  aA  b,„,  b,„,b.  .  All  of  these  quantities  are  dimensionless,  so 

that  when  they  are  used  in  Equations  (38)  the  resulting  values  of  6D2  and  5C2  are 

expressed  in  the  same  units  as  the  degree  variances  kn .  The  maximum  value  of 
n  =  174  was  imposed  by  the  core  storage  available  in  the  CDC  6700  Computer. 
The  results  therefore  reflect  the  influence  of  geoidal  height  wavelengths  greater  than 
about  120  nautical  miles.  The  value  of  n  could  be  increased,  and  the  influence  of 
shorter  wavelengths  represented,  by  modifying  the  an  and  bn  program  to  utilize  the 
disk  storage.  This  in  itself  would  be  expected  to  increase  the  already  long  running 

time  of  the  program  and  it  would  be  increased  to  an  even  greater  extent  as  a  result 

of  the  larger  value  of  n,  since  the  number  of  values  of  the  P™  and  QjJ1  that  must 
be  computed  is  approximately  proportional  to  n2. 


XI.  NUMERICAL  RESULTS 


As  an  example  of  the  computation  of  impact  errors  from  Equations  (38),  the 
degree  variances  kn  (in  meters2)  have  been  computed  from  the  empirical  formula 


e-b(n -  2) 

k  =  a - r 

n  nc(n  -  1  )2 


(43) 


with  the  numerical  values  a  =  6510,  b  =  0.00117,  c  =  1.058.  This  formula  and  the 
values  of  the  parameters  have  been  obtained  by  E.  D.  Ball  as  preliminary  results  of 
a  study  now  in  progress  on  the  combination  of  satellite  and  surface  gravity  data. 

The  following  table  is  a  listing  of  the  results  obtained  for  the  trajectory  having 
a  range  R  =  3000  nautical  miles  and  a  time  of  flight  T  =  1920  sec.  Similar  listings 
have  been  produced  for  all  45  trajectories  used  in  the  study.  To  interpret  the 
listing,  consider  the  last  three  entries  (NAV  AND  IN-FLIGHT)  in  the  last  line 
(n  =  2).  The  value  RMSD  =  >J 6D^  =  122.31  meters  is  the  root  mean  square 
down-range  impact  error  that  would  result  if  the  actual  field  of  the  earth  had 

degree  variances  given  by  formula  (43)  (for  2  <  n  <  174,  the  higher  degree 
variances  being  zero)  but  only  the  normal  field  were  taken  into  account  by 
navigation  and  guidance.  The  value  RMSC  =  5C*  =  89.90  is  the  root  mean  square 

cross-ra  n  ge  i  m  ojet  error  under  the  same  conditions.  The  value 

RMSR  =  JSD2  +6C5  =  151.80  is  the  root  mean  square  radial  impact  error.  All 

values  of  RMSR  in  meters  can  be  converted  approximately  into  circular  probable 
errors  in  yards  by  the  expression  CPE(yd)  =  0.91  RMSR(m). 

Under  IN-FLIGHT  ONLY  the  entries  RMSD  =  244.66,  RMSC  =  126.30, 
RMSR  =  275.34  are  the  errors  that  would  result  if  the  actual  field  were  the  same 
as  before  and  if  navigation  accounted  for  this  field  completely  but  guidance  used 
only  the  normal  field.  Under  ALL  NAV  ONLY,  the  entries  RMSD  =  145.64, 
RMSC  =  49.20*  RMSR  =  153.72  are  the  errors  that  would  result  if  the  actual  field 
were  again  the  same  and  if  guidance  accounted  for  this  field  completely  but 

navigation  used  only  the  normal  field. 

The  navigation  only  errors  are  further  broken  down  in  the  remaining  entries. 
Those  headed  NAV  HEIGHT  ONLY  are  the  impact  errors  which  would  result  if 
navigation  took  into  account  the  vertical  deflections  of  the  actual  field  but  failed  to 
account  for  the  geoidal  heights.  Those  headed  NAV  HORIZ  ONLY  are  the  impact 
errors  which  would  result  if  navigation  took  the  geoidal  heights  of  the  actual  field 
into  account  but  neglected  the  vertical  deflections. 


55 


r 


t 


3  -4  «  -ft  3  n  (KThwian 


^  ft  9  »>•  •9«*«(Mn4»9^9»' 


«4«4^(ViftJMISj«Wft4Ml«ftnlftfftr)mift^«^  r  »  4  ft  -f  »IMMA)A«tf«Ar<MA 


<£  o  3  «*>  M  •  ft  o«\iiAK.'nci^NfwMl^fu*<'3(^«ft.ift  .»  m  o  3  o  •#  u  Q  49  ">  o 

tvi««f  rn  O'  cm 

t  r  ft  .#  ft  4  r  r  r  J  ■»  ft 


0«)")KB9>(r  0  j  4  MJ1  I  H  C  -I  H  «  t  ^  *  ft  »  IP  9  N  •«  O’ 

jjj®«noHM^w  i  •  k.  N.  &  -a  tt  «4  rj  rw  *9 

H  H  H  4  w  «*  *8  «*  -4>4^NNMNNNNMMN<Ni^MNAJNn*)M)1  ft  « 


-«  -4  f  tftJKftJ  3  ^  >J  f  o  tftt  tM^o9s  0  #  ft  ^  o 

K-  O  %i  ft  t0  h.  O'  O  «•«  ^<  >U\<0S«KaKjftn-tlA  0ft«9O^(ylftft4lA  Oft*-’* 

•<w«-^^-4’'i>>irufVjftJMftiftift/ftft»,«ft*ft»ft^l»ft»ft*ft-ft  ft  ft  ft  t  ft  ft  r  4  -f  ft  ft 


*>  O  If  t  3  OftlK.ftJ».««iDn  ft  «  M  0  O  ft  »!  'J  ft  ft*  "ft  8  <3  ft  S  ^  ft  *»  *8  .ft  ft-  *"«  ft 
(M  4  t  V  ^INNO«IT>0>ftOHftH^^1  ^ft  ft  4  r  aMA  i£  <0  O  »  lv  h  «  0  «  (ftft 


*  ~4  -4  M  ft  ft'O'ftfO'ftrO'ft  ft  -4  ft  ft  ft  ft  ft 


r  ft  .»  IT  o  ft  «4«  3  IT  'ft  *4  ft*  |ff4^  |H\JU'«i-^  ft  N.  ®»  «ft  I^)k  J»Hl\)4U'*'v« 
o  ,\  0-4  ft  3  «a  fw  fta>»»ft*'afvjft»r\i>«ftaiNjft>  *  lA  K  o  a«  3  3  "9  ft  to  >©  K  ®  f 


ftWMlMivrv/Nftftftftftl  ft  4  ft  ft  ft  ft  fttMAAIMA^tMT  0  0-0  0  0  ft  ft  ft  <0 


r  <n  -4  a  c  «  (T  ff*  s  N  (T  j  o  ,.M*  «  f  oft*.  f  M  ft  4>  O  a  N  *0  AftKelTOHN 
3  ft  x>  .*  "8  ftKft»<ftliftftOftHN4i(\K0»O(Mft  ft  AK  t  0*^  O 

333f\jf\j*\lf\if\jWK»Wfn-«lft4  4  4  ft  ft  ft  /aMAnftftftAftftOftOftftftft 


^  ^  c  »  -  j  h  ,>  j<  w  »  v4  ro^jr.  .jftftinyjAKoftuiftofta  o  «8  *'*  i*\  «  o 
M  ■**  ft  ft  .8  f\.r  J  Oft'^0*5*'flT'fa033-4-4-4fVir\J<\j\ir'l-»>3*8  ft  4  ft  ft  > 


vWin  ft  a  <8  (A  ft  ft  ")  aK  1  e  t  «  8  K  rtj  K  n  9  ft  «  O  N  ft  D  *>  (T  -4  <\J  "I  ft  Oft 

r  a. i  0-4  r4)ejN40'»J'o(u9tf'Jftir  JV8iftft«ftDHiv  ^  aan-  o^ 

/>  3  -4  3  f\i  JI\,(Vj*n3~13*n-0  f  J  ft  ft  ft  ft  ft  ft  A  1ft  A  .»  \il  .1  j  .1  J  J  Jft  ft  JO 

r 


-3 

4  ■  O  O'  <f  >-||i'U(T'  J  d  IT!  (T  3  3  ft  »  (T  N  ft  o  O  1  "J  18  y\  4  N  O  T  O  H  M 

SJ  J  Ij  nH8,tsH48  IflOOHN  »AN«83N8  ft  A'»UI«3HM8t  JN8 

*  A  ^4-4-4  •ji'IA/*jN88818ftft»ft  ftftfiftlAAAA  l  ft  0  >00003090 

3  r 

T  d 


N  •<  »  t  8<  ft  «T  8  ft  O  r  ft  «4  tft  O  8  Jft  W  ftftlTAjftftftlM  |  ft  {A  H  ft  0  «  >4  8 
M  8  3  ft  rir  ,%  0  0  **•  —  N  f>  O  0  3  -ftl  ^  *  O  ■»333-*33'V"SJ\i!-'a3-0-0  3  ft  ft 


0  ■"  3  O  r  .3s-ma*8(\J*3ft  -\0^O*33f'J'O  ft  ft  ..,  3  ft. 

■»33-4-433.4»\JfM‘>4'\jr\jr\J'\j'V('g\J»8,8’*,)'<>33*833  1  3  ft  ft  ft  ft  ft  ft  ft  ft 


n  8  ■*  3  a  i  a  3  »  o  o  3  a  3  503010  a  333000  1  3  3  ij  03  3  3  3  a 
J  ’  o  o  0-3333  33303330  3  333303  333  333  13  33303  3 

03333  3  03333303303300  3003333300  »  3  3  3333 


O  *  3  ■**  ft  0  Is  ft  IHM^ftAftsafta^NMftA  0S.«<ff33M3  ft  .9  0*. 
n-»HHH4-4M(MMNNMPgNMN88889888889ft  ft  ft  ft  ft  ft  ft  ft 


O  A  ft8M8 


<f  • 


i 

i 


1 

1 

1 


56 


TO*  rj 


I 


-4  «A  <ft  f  L 


IKOffiOHI 


aojaaooo  iOH^^«4>*Mi4HHiyMNNNNiVIVN 

- - - - . 


r«<  hi 

PM>. 

i  o  o  rr*  &,  m  »  p  p  iji  p  j  ^  ffji  u 


<o  ®  t»  n  «o  a 


®  o  o  <n  ©  m  m  i 


t  -•  j*  «c  no  a  o  ow  h<tk  ft  r*j 

-4  0.  \j  K  O  09>e<INnn  HfUM-tP'lHNntl()|S49fli«NfDilsV»MtO 

ft  .O  0  ^  0-0-0  O  a  -0  0  OKKKkKNSKKNf  • 

•4  «4  «4  «4  «4 


fft  t  9«'l*)UN0f44Sol«hBHlH0M44Ha4n4<l9**9aN0K1 

HMK1  j  If  0  *>  *e>-4M«1tr».0k.X>O«4<\j,*tft  ft  *-  (T  -4  N  .ft  ft  k-  *  'V  “ft  ft  k.  Cft  *4  nj  Mf  «  w^iftKO 


mKp(KT9'j'0>iftnooioo^«<H4*«HH4NNVNNNl,l1,»«tt  9  0-#  9  ♦  ft  III  III  ft  0 


4  OMK  H  30  ft  <T>  J  T*  ft  <T  fO'f0>4O'ftHl0fge  tQ0N9lAM9«0iaN|0H«4^k,  B  A  tHH 
i  i  f.  S  fl  «  Jt  n  -AO  j  .4  pg  rg  *T  rft  MiA0dkSff0feoHNA(»)4-A>A0ft0«I<«*<4UK|4 

t  t  0  0  »  *  >0  I  1044  ^444  AJ^Vi  ,MAIA|0ir>0|{l  >1000041  00 


n4i«i*>t>ol«'’'  i  r  H  M  -A  t  9  *•  »  4  A  n«HtknikN0«4IAH09JBIANA^  (ft  ♦  »A  ♦  ill 

>  I  ft  O  >-  o  n  a  j  ->  f  x  j  r>.  IO't-MHji.xs  B»h^A|,i  0»1>  4-g  »  Ak^ft-g  0  o  ®  « »  'M  4  41 

o  o  o  o«i»  o  ft  ft  o'  cr  tr>  < y*  ft  o'  -ft  j  a  »  iDoj!jgH4»4wft4^ivft(^M''jKiw'AiA"i  r  0  r  0 


.3  44  -4  ft.  fft  I  Ml  ft.  o  n  g  l\J  J  ft  A.  0'ON4KAH  M?iNlllBH<BH|(.}'  0  •  "»  0*  4  ®  '*»  (ft  ft  «4  ft.  f 
TOHftJ*ft  r  ft  X>  Cl  ?  -4  N  «  riOft  B0OH*|  J  IA  fl  B  9*  O  M  ft  f  9  •.«e»-4v>4^>ft-  -ft-»rvi0,f\ft- 


HNIMN\M9J  04  9  9  9  IT  ft  ft  l/l  'ft  4A  IA  vf*  9  9  O  >0  0  O  k  n  s  ft  -> 


3  .4  -«  r\J  ig  ift  0  .>4-  0ft  <0*  »  »4  ftf  fill  JO  a  M  Kft  >«nNlftK<9A)09>ft9|O(ftl*>>0CTtn(ft  9  Cft  0  ft  U' 

OJ  .  C'ft^-.fsim  t  A  3  k.  -ft  CJ  -4  fti  *ft  9  ftfteftra^ft^jAfl^tfto-iAO  r  />  k.  «  ^  -4  rvj  9  ft  ft.  w  a  <M 

4H-^rjlMMiMMvyij'\jfti1ftAft«<>')ft43#0|3  09ft  U  ft  M>  ift  Jft  ift  4  if  0  O  0  0-OftN- 

-4  -4  -4  *4  4  4-4-4  4^^4  444444.4,14444  4,44*4444H^4^^l44W»4  -4  -4  -4  -4  4-4  >4 

r  .p  0  U  f  -  o  O  j  O  -4  J  |  ■*•*»  jo  9  0  IT  ft  tfU1  K1  OO  r  «J  H  y,  T  ft  OM  I)  4  ,1  O  ,l  o  IA  -J  O  4  n  n 

a  O  J  .  ft  ft  O  O'  (ft  Ati004444MNM»>*liJ  9  IT  >lft»0  0^ft«t'ft<ft  ft  3  ft  4  M  M  ft 


ift  0  ft»  •  ATOr\|ft  r  ft  Je04Ml0«O9(  f  0«Hft  0  ft  ftj  |ft  «  f\J  |\0  Ak4  A  ft  IP  ft  IP  ft  0  MO 
O  O'  Cl  ^4  (\|  *ft  ft  «.  o  T  a  4  M  J  .10ft  a  U  4  M  ft  f  O  N  ^  T  4  M  ft  ft  «ii  «  ftcft*\jw>.ffti0»«ft(ft44f\|0  ft 


i-4  g  ftj  -\J  (\j  r\j  rg  ftg  - 


ft  AO  ft  rft  -ft  -ft  .»  f  f  >t  9  ft  99. t 


I  Ift  •>  ft  >jl  .ft  J  0  0  O  O  0  0  '  -  0  ft  ' 


j  4  4  m  M  ft  0  .Ift  0  0ft  ftaHftj  4  ft  i  O  a\f  ft  ft©  nkaftJTMy'Tft  *>  o  ift  (7*  »  t  f  J>  'P  o 

4  .  O  ^  -14  M  -ft  f  ,0  k  T  J  4  PJ  ft  pl'IWft  f  m  k  O  O  4  AJ  0  ft  k  OOftJ 


(4H4M^A|MMN4  \J  \i  ft~9>ft  t  0  9  9  9  0  0  9  Ift  ft  'ft  (ft  Ift  ftftOiOiA  O  0  9  9ft-* 


■ft  P  flo  N  t  ft  ft*  -4  >ft  O  "ft  3  "ft  ft  k.  ft  4M  ft  k.  3MlPkOM1p«a’9  0«4^kOft  9  «  4  Jft  0  *4  0  0«.  «4 

4  0  #  lit  Ift  IT  O  0  ft  09  p.k.9.  O  O  0®AA4ITftann4444(Mr'JMft"9,ftl*  p  f  0  'ft  ft  ft  0 

M  *\»  fti  M  g-vj-gft.  ftiftiftj-viftjftjftjrg  gfti  vjPj'MNMMftftft'ft  ft»n,ft>ft-ft'ft'ft»ft'ft,ft’»»’0'ft*ft'ft  ft  ft  -o 


M  ft  '  ft  ftftj  r  J  O  J  "\J  »  J  O  O  ftj.ft  -  IMjIkOft  O'AVI  ft*M  AI0»I  -  4  iftl"  4  »  ft  ft  ft  m  f  -» 

4Tinn44444-\JMMMM,9ft-9ftJ40  tiftft  AiP00  0f«kkefl94iAoO449fAift  ♦ 


_4_4-4-4-4-4-4-4-4  4f4H44444H4444444444444H44HH4HN  NJ  U  U  *U  (\J  IU  >J 


ft  i  »  r  i  i  ft  »  i  ft  ft  >n  to  jaoo  1  l  itftffto^ft  jo  inoaonftao  »  ft  ft  ft  n  -j  o  irj-ft 
ft  ift  i  ft  9  n  ft  iftftcjuftftoftoftaft'iftftftftoftftoftftft'iftaoftftricj'ftaftftft 


J=3fft  ft  ft  -  rl  .ft  n  ft  ft  ftftlOl  lftOOOOAftOOlftOOBflftOftOOOftO'laftlOft 


vJ  -ft  OftftJ  9  30  ft  ftj  0  o  ©  •a  ftj  ■>  k-  ft  NJ  .ftk.'ft’H  o  ■*  ftj  ft  0  M  lA  »  A  k  *4  ft  4*  o  ft  o’**  A  »"» 

ft  ft  3'-»'3-4-4-4-4-4f\J'\Jftjr\J\J*ft*ft  ft  ft  0  0  M  PUI  n  A  »  0  SkkMI  BT'0Aq{»44'g'M,i  0 

4-4-*  4  H  4  4  <4  4  4  4-4  4444444444H444444444444444VMNM\MMM 


-4  ft  ^  Ok.  Oft  *-ft«g.4cftft«t>9.  On  0*AM-4a!?'«N  0  ift  f<IM4aT0k  o  ft  9  "ft  kg  -4  Iff  «k  0 

«AP9009nrB90.'.  .k.k  k‘,  1.3303999999  lift  ft  ft  ft  >  ft  ft  i  f  0*  * 


i 

■i 

58 


I 


*  99  9*  Q 

♦  O  A.  O' 


O'  «  S  f  H  f  O  D 

t  O  e  o  N  }  >0  9  H  n  .a 


•  w*  9  k  « 


4  a  «  IPK  « 
•  t  n  ro  if)  >4 


n  n  un  m  .0  o  r  o«NKKN.KC«««99f?iiS<«^<4MM«l44iniOK*oN4«^li' 


•  iP>-  >  k  *  *»  « 

-i  *4  u  9  o  .f  m  nAN«e'T>io^i\«n«i4wwtf»>(>.C0'9e^«4Ajn 


-4  .4  -4  t  -*  m  ®  CM  3  0  <** 

.4  -4  -* 


O' 

„ 

~t 

p 

j>  M 

-4  ~S 

O'O'OCMil  O  9 

-4  m 

K  44 

44 

*  JJ  JJ  ° 

■4  N 

44 

9k 

9* 

*■* 

Ai 

P  A 

O  o 

-« 

o> 

i  »  1  ^  «N 

o 

*9  DONilM'.NJ'gMlO* 

O  A 

O  « 

1# 

f  lO  J  «N9 

3 

, 

•4 

*4  -4 

*4  (SJ  fVJ 

IM 

M  fg  n  ri  9  *) 

♦ 

i  (4«iAIAM».N««*ai3 

*4  PM 

900 

40  O  9  9. 

ii>  m 

•M 

.4  -4 

•«•  -4 

"* 

*4  <4  H  *4  M  -4 

** 

ij  PM 

N  CM  IM 

M  99  >9  9 

t  -£>  <M 

-4 

O'  9. 

43 

'40ss9IAi\t(U 

9  «- 

O  0 

a 

Id  4  «  t 

0  0 

t 

9J 

p 

**• 

O 

CM  ft 

®  »4 

t 

tr 

4  IA  C  Nit  4 

•A 

o  a  i  jN9'hNo*«oh  r  « 

m  n 

«  K 

M 

in  IM  44  K 

9 

o 

9 

vC  N 

•*.  9 

Ik 

9 

« 

O'  O'  O'  O  O  44 

N  N  O  9  •»  *»l»  0NN«»  HN« 

n  K 

(P  ^ 

0 

»  0  «»  «o 

O'  ^ 

•ft 

-'4>4f0MNcgMVyMN-MNNNNNN',>*’M99 

99  JP 

3 

J\  ft  ifi  9. 

O'  P  K 

** 

IM 

9 

P 

ir 

*r> 

O 

iKffWff'K.N*  iN  J  i4  9  N  3 

o  a 

0*  99  M 

44  99  (M  cm 

n  «n 

»3 

» 

D 

p 

« 

<r  ci 

rv  9 

t 

o  ^  «Q  NltA 

»4  ^  4tA(r  /  (MAN 

0*  O' 

O'  m 

CM  M  0  P 

K  4 

9 

J 

o 

0 

j  f- 

r 

f 

••«  ►.  ®  ®  •  ® 

« 

fl  Pfff000044  4  CM  CM  O*  i 

J  A 

kD  « 

3 

CM  k.  1»  vb  K  IT 

O 

*4  >4 

•4  -4 

CM 

IM  CM  IM  9 

J  0 

CM 

, 

M  a 

a  r 

M  f  X>  Pn  3  O 

in  n 

9  K 

44  *9  O)  m 

•9  -» 

•n 

D 

-* 

i 

O  3 

m  D 

M 

1  O  -4  ,1  l»  CM 

4> 

*4, -to  \  0  0M  n^>  1  9  N  9  p  , 

<M  <T* 

0  in 

0 

IN  P  1 

O  9 

0 

J 

ft 

, 

xO  ■£> 

0 

*  r>  q  o  cr 

CT 

rjn44Nlj99p  0  0  K  «  O'  D  ^  9 

0  0 

n  0  b  k 

•O  K 

p 

-4 

•4 

4  4  r«  4  r<  H 

■4 

M  IM  'JMf'JMCVJMCM'NJfMfM-MM>9’OK> 

99  19 

t 

t  m  n  0 

O  M 

t 

CM 

O' 

r> 

9  m 

O) 

f  xC  O'  M  4)  -4 

P  X4NJB9444aOiAU'NN  P 

0  *- 

9- 

O  J  *9  O 

O  *9 

CM 

T 

M 

p 

0  « 

rs  CM 

r 

11  n  N  ,3  K  O 

'j 

•  P«\i)4iOHS  P 

44  3 

CO  IM 

'» 

0  M  0*  0 

9*  P 

► 

fci 

■9 

9 

9  ®  O'  O' 

O'  O' 

oii  i  n  o  h 

•4HNMIM919  P  -p  in  n^>  OK 

m>  it 

O  »4  CM 

f  K  O  t  it  V 

9 

*"* 

H 

“* 

** 

4  -4 

•4  -4 

-4 

-4  "Vi  <M  \J  M  %i  CM 

NMCMNN'JMMNNNNNNNNN 

9  <**  *9 

•99ft 

ft  K 

19 

>1 

9 

O' 

a 

P  4 

»*  o 

O' 

A  U  ,.  O  M?  4 

J 

«  n  s  j  >4  <r  ro'CMUJ-iO'crcM® 

•  IM 

44  K 

IM 

ao  0  A  K 

SO  -4 

O 

1 

0 

o 

3  cm  -»  m  p 

O' 

•4  9.A',  r  CM 

■t 

J  ip  4  4  N  *  CM  01  O'  CM  0  o*  K»  O  M  K  9 

91  m  *9 

9N  J  1 

C9  IM 

CM 

9- 

„ 

►k 

K 

®  *3 

o  ® 

o 

e 

tr  o*  o'  o*  <r  * 

a  J  H4  4  4Pi\rj"»"19  f  i  1/1  A  J 

O  K 

« 

cr  0  -4  i9 

•ft  O' 

O' 

•4 

“* 

** 

**  -4 

44  »4 

’■* 

44H444(\jNNN\Nc\<N^NN,VNNNNNi\jN 

<M  M 

CM 

fti  *9  •«  9 

9  9 

f 

n 

’9 

«T 

IT  -4 

c 

^  er  o'  4i  o  vO 

ll  -4«5r.fk0'4Ooj44in-4  .  A3 

«  n 

j»  a> 

O 

m  j\  .1  p 

fM  K 

p 

■ 

•> 

i 

l>  p 

Cl 

■4  (\i  <»1  (,*  If  k 

<r 

ci  cm  cn  ir  K  O*  (M  n  O  4  ,p  o  p  O'  o 

co  fM 

9  K 

,  •-  4.  » 

19  CM 

J 

| 

i 

r 

, 

P  P 

-1  » 

, 

L 

l.'  .>  »n  it  >  to 

flj)fljjiOSKN«t)«j(I'IJ'3 

44  CM 

C9  f 

03  -4  ■£  CM 

■i  0 

"* 

44  -* 

•4  44 

44  CM  M  9 

r  J 

» 

P 

T 

9  *4 

4»  «1  K 

►- 

^  k.  «  o  <M  n 

•o 

"9  «coeeKK.(T(M  O  C9  44  »M  -O  C9 

»  cr 

ca  n. 

t 

CM  O  P  O' 

♦9  O 

K. 

* 

" 

a 

M 

3  J  K  O' 

rrx 

i\  n.  «  ^  *  3 

o 

HWJirHjKO  t  K  «4  41  "  -9  o 

19  t> 

in  4 

n 

®  0-41- 

-4  a 

IM 

o 

O 

u  T' 

•J  O 

41  O' 

ff>  T»  **  3  3  3 

HHH4  JCMCM-9  9  9  P  P  J  n  > 

■D  xD 

-  w 

•0 

1*  O  M  9 

b  O 

O 

”* 

4 

** 

*“  *" 

*4  44 

-4 

OhoNNNNNNNNNNNNNNNNNNNNN 

IM  CM  9J 

CM  9  *9  99 

9  » 

1ft 

v) 

v, 

o' 

t  4 

O'  o 

ji  n  0  s  A  4 

O  cm  K  P40'3'J»Nii)H(riJ*N  o 

W  CM 

44  K 

3>  O  ft  ■» 

43  -o 

■1 

9 

' 

** 

i) 

.3  M 

"»  .A 

*- 

n 

-4  4  ,1  K,  »  ^ 

f 

0  PH  PC^O*>j  H-IJ  J  P9  ON 

N.  -9 

■9 

44  CM  ^  O 

9  "J 

CM 

p- 

„ 

•a 

o  o 

O  41 

n 

•» 

«  »  *  o  *n 

■a 

■»-*H-4H44*vi-\iAJ«N"4  9  9  P'Ain  0 

b  . 

0 

04  IO 

ft  O' 

O* 

-* 

-4 

*4 

4 

-4  -4 

-*  *4 

-4 

-4 

44H-4HN'\iNM'VIN^NNMAiPgNNNMN 

■M  Ai 

M  M 

91 

•M  9  9  *9 

19  9 

P 

P  N  -4  I  O  CM  L  cr  ->  K  H  0O  fJi»1M.l'J(T'iArtK  tlh  AlriOUON 

1  1  f*  K  ».  O  O  O  IT  *>  n  jww4MNfO^^  P  ft  3  OK®«''’l'*Hl\i"l  #  ft  QK.'T^  9  ift  i1  .f  iM  9 


i-no-o  990999  I  ,  r  f  p  ♦  f  *  »  -t  f  -f  t  -t  J  -t  t  t  MTiA  >  U\  .n«  1  A  On  0  J  ^  ft  a 


,i  *o  9  »v  pftj.-*0'nr>T.»'-490*  0  \|  »«*.<  ^  4  y  -3  »  i«  n  K 

f  ft  o  O'-  a  tr>  tn  ^  m  -I  ft  .0  k  9  iMlirsai\hA*>4u.  o  ")  •<  o  i*  (\j 

M-y"'ji'J,vJfy'\jM  9  9  o  -<1  r<j  1  1  p  t  p  ft  ft  ft  ft  'ft  .0  3  a^h»ort'c3aN"3j\^a  i  fv  j\ 

-4-4  ^  H  -*  >J  I\J  *0  P  1  t 


win  »  i  rj  i  ci  a  i  n  i  b  a  ^  ji'Jarjmdrj'ioioOTO  isobtiii  ■»  a 
m  i  n  -s  ^  i  n  ^  iTnnai'no-jain^cJOo^'ino  i  i  rj  n  i  i  -t  t  n  -s 


3  1  in  ^  o  o  n  i  u  n  liana  jaano-jinnaann  moan  s  n  n  n  a  m  a 


59 


RMS  RADIAL  IMPACT  ERROR  IN  METERS  DUE  TO 
NEGLECT  OF  ANOMALOUS  GRAVITY  IN  GUIDANCE 
AND  NAVIGATION 


FIGURE  6 


FIGURE  7 


RMS  RADIAL  IMPACT  ERROR  IN  METERS  DUE  TO 


FIGURE  8 


For  the  preceding  rows  of  the  listing  the  actual  field  is  again  assumed  to  be 

that  represented  by  the  kn  for  2  <  n  <  174  and  the  impact  errors  result  from  the 

omission  in  the  gravity  computations  of  terms  of  degree  n  greater  than  or  equal  to 
that  indicated  in  the  first  column;  this  omission  occurring  everywhere,  in  guidance 
only,  etc. 

The  radial  impact  errors  for  all  45  trajectories  resulting  from  the  omission  of 
all  terms  n  >  2  are  represented  in  the  error  contour  charts  of  Figures  6-10.  For 
given  values  of  R  and  time  of  flight  T  the  total  radial  impact  error  (in  meters)  and 
the  contributions  to  it  can  be  read  from  these  charts. 

An  interesting  conclusion  concerning  navigation  errors  can  be  drawn  from 

Figures  9  and  10.  The  sketch  of  Figure  11  represents  the  superimposed  40  meter 
error  contours  from  these  two  charts.  For  points  to  the  right  of  the  more  vertical 
contour,  the  neglect  of  geoidal  heights  by  navigation  produces  radial  impact  errors 

greater  than  40  meters.  For  points  to  the  left  of  the  looped  contour,  the  neglect  of 
vertical  deflections  by  navigation  produces  radial  impact  errors  less  than  40  meters. 
Hence,  in  the  shaded  area  the  neglect  of  geoidal  heights  makes  a  greater 
contribution  to  the  impact  error  than  does  the  neglect  of  vertical  deflections.  By 
plotting  the  intersections  of  corresponding  contours  in  Figs.  9  and  10  and  joining 
the  points  with  smooth  curves,  as  in  Figure  12.  a  region  is  delineated  within  which 
the  neglect  of  geoidal  heights  makes  a  greater  contribution  to  the  impact  error  than 
does  the  neglect  of  vertical  deflections.  This  region  covers  a  large  part  of  the  area 
representing  possible  combinations  of  range  and  time  of  flight. 


R  (Mint  Mt  ) 


FIGURE  11 
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At  present  the  navigation  system  is  compensated  for  vertical  deflections  in  areas 
where  the  data  are  available,  but  transmits  to  fire  control  no  output  of  height 
above  the  reference  ellipsoid,  for  which  data  on  geoidal  heights  would  be  necessary. 
Fire  control  makes  some  allowance  for  geoidal  height  in  the  computation  of  target 
offsets  due  to  tesseral  gravity,  but  makes  no  allowance  for  variations  in  launch 
depth. 

The  height  of  the  submarine  above  the  ellipsoid  (or  depth  below  it  if 
height  <  0)  is  just  as  truly  a  navigation  quantity  as  arc  the  coordinates  of  latitude 
and  longitude  and  should  properly  be  accounted  for  by  fire  control  and  missile 
guidance  as  one  of  the  three  components  of  initial  position,  instead  of  being  only 
partially  accounted  for  in  the  target  offsets. 


FIGURE  12 


.1 


1 


It  would  appear  that  a  comparatively  minor  modification  of  the  navigation 
system  would  allow  it  to  produce  and  transmit  to  fire  control  the  height  of  the 
submarine  above  the  reference  ellipsoid.  Figure  13  illustrates  the  computations 
involved  in  a  three-channel  inertial  navigator  which  employs  data  from  the 
electromagnetic  log  to  control  errors  in  the  horizontal  channels  and  data  from  the 
depth  gauge  to  control  errors  in  the  vertical  channel,  the  latter  in  such  a  manner  as 
to  produce  an  output  of  height  above  the  reference  ellipsoid.  The  method  shown  of 
introducing  log  data  is  identical  with  that  now  employed  in  the  MK  2  MOD  6 
SINS.  Also,  the  method  now  in  use  to  introduce  vertical  deflection  data  into  the 
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FIGURE  13 


SINS  is  practically  equivalent  to  that  indicated  in  Figure  13,  the  stored  and 

interpolated  values  of  the  deflection  components  being  equivalent  to  the  gK  and  g^ 
outputs  of  the  Gravity  Computer.  The  MK  2  MOD  6  SINS  also  includes  an 
incompletely  instrumented  vertical  channel  which  is  stabilized  by  the  depth  gauge 
and  produces  an  output  of  vertical  velocity  vh  which  is  transmitted  to  fire  control, 
but  no  output  of  height  h. 

This  incomplete  vertical  channel  can  be  regarded  as  an  approximation  to  that 
shown  in  Figure  13.  In  the  approximate  instrumentation  the  Geoidal  Height 

Computer  is  omitted  (equivalent  to  taking  N  =  0)  and  the  output  of  the  second 
integrator  in  the  vertical  channel  represents  approximately  the  height  of  the 

submarine  above  the  geoid,  but  is  not  transmitted  to  fire  control.  By  including  the 
Geoidal  Height  Computer  in  the  manner  indicated  in  Figure  13,  the  output  of  the 
second  integrator  would  represent  height  above  the  ellipsoid  and  could  be 

transmitted  to  fire  control.  Use  of  this  value  in  the  initial  position  computations 
would  account  in  the  theoretically  correct  manner  both  for  geoidal  heights  and  for 
variations  in  launch  depth. 

The  Geoidal  Height  Computer  would  in  practice  consist  merely  of  a  means  of 
storing  and  interpolating  geoidal  heights  at  a  set  of  grid  points,  exactly  similar  to 
the  method  now  used  for  the  vertical  deflection  components,  and  could  be 
implemented  by  storing  the  geoidal  heights  on  the  same  tapes  as  the  vertical 
deflections  and  subjecting  them  to  the  same  interpolation  process.  Furthermore,  the 
basic  data  required  to  produce  the  table  of  geoidal  heights  is  identical  with  that 
required  to  produce  the  tables  of  vertical  deflection  components,  and  the 
computations  are  very  nearly  the  same. 


XII.  SPECIALIZATION  TO  A  PARTICULAR  LAUNCH  AREA 


In  introducing  the  statistical  treatment  adopted  in  this  report  it  was  admitted 
that  a  more  useful  statistical  description  of  the  impact  errors  would  be  one  for 
which  a  single  target  is  considered  and  the  launch  points  are  uniformly  distributed 
over  some  limited  launch  area.  It  seems  possible  to  utilize  the  statistical  treatment 
which  has  actually  been  employed  to  obtain  results  which  are  in  some  respects 
representative  of  what  would  be  expected  in  this  more  realistic  situation. 

This  extension  rests,  first,  on  the  observation  that  the  initial  position  errors 
depend  on  the  geoidal  height  and  the  vertical  deflection  components  at  the  launch 

point  and,  second,  on  the  fact  that  the  in-flight  errors  depend  principally  on  the 

geoidal  heights  over  an  area  which  does  not  extend  very  far  from  the  launch  point. 
Three  reinforcing  factors  contribute  to  the  truth  of  the  latter  statement.  First,  the 
effects  of  the  short-wavelength  variations  in  geoidal  height  which  might  distinguish 
one  launch  area  from  another  are  attenuated  more  rapidly  with  height  than  are  the 
effects  of  the  long-wavelength  variations  and,  hence,  would  be  most  significant  near 
the  launch  point  and  near  the  target.  Second,  disturbances  of  the  gravity  Field  near 

the  launch  point  make  a  greater  contribution  to  the  impact  error  than  do  similar 

disturbances  near  the  target  because  the  missile  is  moving  more  slowly  near  the 
launch  point  and  is  affected  by  the  disturbances  for  a  longer  time.  Finally,  a 
velocity  error  produced  by  a  disturbance  early  in  flight  has  a  longer  time  to  act 
than  does  a  similar  velocity  error  produced  late  in  the  flight. 

For  these  reasons  it  is  plausible  to  assume  that  the  impact  errors  would  not  be 
much  different  if  the  disturbances  of  the  gravity  Field  were  limited  to  an  area  near 
the  launch  point  than  they  would  be  if  similar  disturbances  extended  over  the  entire 
earth.  Hence,  the  statistical  treatment  which  has  been  employed  in  this  report  can 
reasonably  be  interpreted  as  applying  to  a  specific  launch  area  by  employing  a  set 
of  degree  variances  kn  which  would  characterize  the  geoidal  heights  over  the  entire 
earth  if  the  entire  earth  were  as  “disturbed”  as  is  the  launch  area.  There  are  several 
objections  to  this  procedure,  some  of  which  are  unavoidable.  First,  the 
short-wavelength  variations  in  geoidal  height  would  produce  navigation  ve.ocity  errors 
which  might  be  more  important  than  the  navigation  position  errors.  Although 
navigation  velocity  errors  have  not  been  considered  here,  it  seems  possible  that  they 
coulj  be  included  and  this  objection  avoided.  A  second  and  unavoidable  objection 
to  this  method  lies  in  the  fact  that  the  statistical  treatment  employed  in  this  report 
necessarily  leads  to  the  conclusion  that  the  mean  impact  error  is  zero,  and  this  is 
almost  certainly  not  the  case  for  missiles  launched  toward  a  single  target  from 
points  of  a  limited  launch  area.  Finally,  and  also  unavoidably,  is  the  difficulty  of 
quantifying  the  idea  of  “near  the  launch  point”  and  obtaining  a  set  of  degree 
variances  which  would  adequately  describe  the  geoidal  heights  in  such  an  indefinite 
area. 
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It  seems  worthwhile  to  digress  here  to  emphasize  some  fundamental  differences 
between  the  statistical  methods  used  in  this  report  and  the  more  elaborate 
techniques  which  have  been  used  in  some  previous  studies  of  missile  impact  errors. 
These  other  studies  have  extended  to  the  treatment  of  in-flight  errors  some  methods 
which  have  been  widely  employed  in  navigation  error  analysis.  In  these  methods,  the 
incomplete  or  approximate  knowledge  of  the  geoidal  height  function  (or  other 
gravity-dependent  quantities:  gravity  anomalies,  deflections,  etc.)  has  been  represented 
by  modeling  the  geoidal  heights  (e.g.),  as  a  stochastic  process  on  the  sphere  or, 
more  often,  in  the  plane.  This  view  is  expressed  in  Heiskanen  and  Moritz’s  Physical 
Geodesy  by  the  statement,  “It  should  be  mentioned  that  the  mathematics  behind 
the  statistical  description  of  the  gravity  anomalies  is  the  theory  of  stochastic 
processes.  The  gravity  anomaly  field  is  treated  as  a  stationary  stochastic  process  on 
a  sphere  •••”.  More  recently  Moritz  ( The  Role  of  Statistical  Techniques  in  the 
Determination  of  the  Earth’s  Gravitational  Field,  International  Symposium,  26-30 
November,  1973)  appears  to  have  modified  this  point  of  view  by  saying  “Therefore 
Moritz  proposed  to  use,  instead  of  an  interpretation  as  a  stochastic  process,  an 
interpretation  in  terms  of  a  covariance  analysis  of  individual  functions  •••  .  This  also 
takes  into  account  that  in  reality  there  is  only  one  individual  earth’s  gravity  field 
and  not  a  phase  space  of  many  such  fields.” 

The  elementary  statistical  methods  of  this  report  seem  to  be  more  consistent 
with  Moritz’s  recent  views  than  do  those  of  many  previous  studies.  In  particular,  the 
degree  variances  are  regarded  here  not  as  statistical  quantities  (ensemble  averages)  but 
merely  as  a  partial  description  of  a  completely  deterministic  (although  possibly 
unknown)  geoidal  height  function.  The  only  properly  statistical  concepts  employed 
here  arise  from  the  consideration  of  an  ensemble  of  launch  points  and  targets. 

Thus,  it  seems  impossible  to  obtain  a  completely  satisfactory  statistical 
description  of  the  impact  errors  occurring  when  missiles  are  launched  toward  a 
specific  target  from  a  limited  launch  area,  except  in  the  case  when  the  gravity  field 
is  completely  known  and  the  statistical  description  is  unnecessary  (cf.  Heiskanen  and 
Moritz,  in  a  somewhat  different  connection,  “This  we  obviously  do  not  know;  and 
if  we  knew  it,  then  the  covariance  function  would  have  lost  most  of  its  significance, 
because  then  we  could  solve  our  problems  rigorously  without  needing  statistics.”) 
With  the  understanding  that  it  is  practically  impossible  to  solve  in  any  precise 
manner  the  statistical  problem  posed  by  a  limited  launch  area,  it  nevertheless  seems 
possible  to  obtain  plausible  results  by  adopting  a  set  of  degree  variances  which 
would  constitute  a  partial  description  of  the  geoidal  heights  of  the  entire  earth  if 
the  entire  earth  were  “as  rough”,  in  some  intuitive  and  imprecise  sense,  as  is  the 
launch  area  under  consideration.  It  is  planned  to  make  such  computations  by 
modifying  the  values  of  the  parameters  involved  in  Equation  (43)  after  some  of  the 
topics  discussed  in  the  next  section  have  been  investigated. 
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XIII.  LIMITATIONS  AND  POSSIBLE  EXTENSIONS 


The  preceding  analysis  and  the  numerical  results  have  been  based  on  the 

assumption  that  the  ship  remains  stationary  at  the  launch  point,  with  the  result  that 
the  navigation  velocity  errors  are  zero.  Retaining  the  maximum  value  of  n  =  174 

that  has  been  used,  it  would  probably  not  be  difficult  to  extend  the  analysis  and 

results  to  include  velocity  errors  for  ship’s  speed  below  about  10  knots.  This  follows 
from  the  fact  that  for  n  <  174  the  shortest  wavelength  of  geoidal  height  which  is 
taken  into  account  is  about  120  nautical  miles.  For  a  ship’s  speed  of  10  knots,  this 
distance  would  be  traversed  in  12  hours  and  the  vertical  deflections  affecting  the 
navigation  system  would  have  this  period,  or  a  frequency  of  about  0.5  rad/hr.  This 
is  so  far  below  the  Schuler  frequency  of  about  4.5  rad/hr  that  the  vertical 
indication  established  by  the  navigation  system  is  essentially  coincident  with  the 

plumb  line.  Hence,  the  navigation  position  errors  could  be  computed  as  before. 

The  navigation  velocity  errors  in  this  case  would  be  simply  related  to  the 
changing  direction  of  the  plumb  line  relative  to  the  geodetic  vertical,  as  influenced 
by  the  speed  of  the  ship,  and  might  be  accounted  for  rather  easily.  The  only 
complication  which  would  be  expected  to  be  encountered  in  this  process  would 
become  evident  in  the  situation  first  considered:  the  nominal  launch  position  being 
the  north  pole  of  the  reference  sphere  and  the  target  in  the  Meridian  of  Greenwich. 
To  discuss  this  complication  it  is  assumed  that  some  constant  ship’s  speed,  say 
10  knots,  will  be  adopted  throughout  the  investigation  and  that,  if  necessary,  all  of 
the  results  will  be  recomputed  for  some  different  ship’s  speed.  Thus,  the  speed  of 
the  ship  is  not  regarded  as  a  parameter  in  what  follows.  An  additional  parameter 
does,  however,  enter  the  problem,  for  as  the  ship  crosses  the  north  pole  it  may  be 
moving  south  along  any  meridian;  that  is,  the  velocity  of  the  ship  may  form  an 
arbitrary  angle  6  with  the  direction  to  the  target,  and  the  navigation  velocity  error 
would  depend  on  this  angle. 

Once  the  initial  velocity  error  has  been  determined,  there  is  no  difficulty  in 
computing  the  resulting  impact  error,  since  this  would  depend  on  the  elements  of 
the  last  three  rows  of  the  matrix  X(T,  t0).  These  elements  were  actually  computed 
in  deriving  the  results  which  have  already  been  presented,  but  were  not  used  to 
obtain  those  results. 

When  an  arbitrary  launch  position  and  target  azimuth  are  considered,  it 
becomes  necessary  also  to  consider  an  arbitrary  value  of  5.  The  complex  impact 
error  p(afi,y)  of  Equation  (36)  would  then  become  p(aJ3,y,8).  To  genera'ize  the 
earlier  statistical  treatment  it  would  now  be  necessary  to  form  ensemble  averages 


71 


over  all  values  of  the  four  parameters  instead  of  over  values  of  only  the 

three  parameters  involved  previously.  It  has  not  yet  been  determined  whether  it  will 
he  possible  to  carry  out  these  operations. 

Assuming  that  the  extension  of  the  analysis  just  outlined  can  actually  be 
made,  it  would  still  not  be  adequate  to  treat  the  case  where  the  speed  of  the  ship 
is  much  in  excess  of  10  knots,  or  even  the  10-knot  case  if  values  of  n  much  larger 
than  174  are  to  be  considered.  In  both  of  these  cases  the  frequency  components  of 
the  vertical  deflections  would  lie  so  close  to  the  Schuler  frequency  that  it  would  be 
necessary  to  account  for  the  transient  response  of  the  navigation  system  to  vertical 
deflections.  It  is  expected  that  the  height  and  vertical  velocity  errors  would  not  be 
significantly  affected  by  the  speed  of  the  ship  as  a  result  of  the  strong  control  over 
the  errors  of  the  vertical  channel  exerted  by  the  depth  gauge. 

It  has  already  been  stated  that  increasing  the  maximum  value  of  n  beyond 
1 74  woidd  be  expected  to  increase  the  computing  time  very  substantially.  This 
increase  comes  about  principally  in  computing  a4n  and  b4n  of  Equations  (42),  to  a 
lesser  extent  in  computing  a3n  and  b3n,  and  not  at  all  in  computing  a2n>  b2n,  and 
a,n.  An  examination  of  a4n  and  b4n  as  functions  of  n  <  174  for  one  trajectory 
shows  that  they  increase  almost  linearly  with  n.  It  therefore,  seems  possible  that 
these  quantities  could  merely  be  extrapolated  to  much  larger  values  of  n  with 
checks  of  the  extrapolation  being  made  for  only  a  few  of  the  45  trajectories  by 
actually  computing  a4n  and  b4n.  The  behavior  of  a3n  and  b3n  appears  to  be  more 
complicated  but,  if  an  extrapolation  cannot  be  employed  for  these  quantities,  it 
would  not  be  difficult  to  compute  them  at  length,  since  the  time  required  is 
approximately  proportional  to  n  rather  than  n2  as  for  a4n  and  b4n.  It  might  be 
mentioned  that  the  accuracy  of  such  extrapolations  could  be  rather  low  because  the 
degree  variances  k()  of  any  reasonable  geoidal  height  function  decrease  so  rapidly 
with  increasing  n. 


APPENDIX  A 


DISTRIBUTION 


Defense  Documentation  Center 
Cameron  Station 

Alexandria,  Virginia  22314  (2) 

Strategic  Systems  Project  Office 
Department  of  the  Navy 

Washington,  D.  C.  20360  (2) 

Attn:  D.  Gold  NSP-230 
A.  King  NSP-240 

Naval  Oceanographic  Office 

Washington,  D.  C.  20390  (2) 

Attn:  T.  Davis 
J.  Hankins 

Defense  Mapping  Agency  Aerospace  Center 
Second  and  Arsenal  Streets 
St.  Louis,  Missouri  63118 
Attn:  M.  Schultz 

Naval  Ordnance  Laboratory 
White  Oak 

Silver  Spring,  Maryland  20910 
Attn:  W.  Coffman 

The  Charles  Stark  Draper  Laboratory 
68  Albany  Street 

Cambridge,  Massachusetts  02139  (2) 

Attn:  H.  Feltquate 
C.  Hrbek 

Singer-Kearfott,  Inc. 

Plant  10,  Dept.  5750 
150  Totowa  Road 
Wayne,  New  Jersey  07470 
Attn:  G.  Blauth 


A-l 


General  Electric  Ordnance  Systems 
100  Plastics  Avenue 
Pittsfield,  Massachusetts  01201 
Attn:  I.  Bender 

Lockheed  Missiles  and  Space  Company 
Department  8115,  Building  154 
P.O.  Box  504 

Sunnyvale,  California  94088 
Attn':  T.  Reynolds 

The  Analytical  Sciences  Corporation 
6  Jacob  Way 

Reading,  Massachusetts  01867 
Attn:  S.  Jordan 

Sperry  Systems  Management  Division 

Great  Neck,  New  York  11020  (3) 

Attn:  N.  Simonelli 
N.  Zabb 
H.  Strell 

Autonctics  Division  Rockwell  International 
Anaheim,  California  92802 
Attn:  J.  Boltinghouse 


Local: 

KA 

KL 

KG 

KG  (T.  Alexander) 

KG  (E.  Ball) 

KG  (W.  Elsaesscr) 

KG  (O.  Schultz)  (10) 

KGD  (C.  Duke) 

KGD  (S.  Young) 

KGR  (R.  Carson) 

KGR  (S.  Gripshover) 

KGR  (J.  Robinson) 

KGR  (J.  Brown) 


A-2 


KGR  (R.  Coppola) 

KGR  (G.  Sitzman) 

KGR  (S.  Lipscomb) 

KGR  (D.  Lando) 

KR 

KAG  (B.  Zondek) 

KAO  (P.  Ugincius) 

KAS  (P.  Young) 

MIP-3  (5)  (5) 

MIP-2 

MIL  (1)  (2) 

M1M 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whtn  Data  Entormd) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER  |2.  GOVT  ACCESSION  NO.I  3.  RECIPIENT'S  CATALOG  NUMBER 


TR-3127 


4  TITLE  (ond  Subtitle) 

MISSILE  IMPACT  ERRORS  DUE  TO  NAVIGATION  AND 
IN-FLIGHT  EFFECTS  OF  ANOMALOUS  GRAVITY 


5.  TYPE  OF  REPORT  &  PERIOD  COVERED 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7  AUTHORf»J 

Oscar  T.  Schultz 


8.  CONTRACT  OR  GRANT  NUMBERfaJ 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Weapons  Laboratnvv 
Dahlgren,  Va.  22448 


1  I.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


12.  REPORT  DATE 

Sept  1974 


II.  NUMBER  OF  PAGES 


4  MONITORING  AGENCY  NAME  ft  ADDRESS///  dill  front  /ram  Controlling  Olll co)  IS.  SECURITY  CLASS,  (of  thlt  topott) 

UNCLASSIFIED 


<Sa.  DECL  ASSIFICATIOn/ooWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (ol  tbit  Rtporl) 

^i^o^1°n,1uimited  t0  u-8*  government  agencies  only;  Test  and  Evaluation; 
jm|^  1974.  Other  requests  for  this  document  must  be  referred  to 
Commander,  Naval  Weapons  Laboratory,  Dahlgren,  Va.  22448 


17.  DISTRIBUTION  STATEMENT  (ol  II la  tbflrtcl  ontorod  In  Block  10,  II  dlttoronl  horn  Roport) 


19.  KEY  WORDS  ( Conttnu •  on  tovotoo  oldo  It  no comoory  «nd  Identify  by  block  number) 


20.  ABSTRACT  (Contlnuo  on  rovoroo  at  dm  It  nocoooofy  md  Identity  by  block  mmbor) 

™is  rePort  presents  an  analysis  of  the  impact  errors  of  a  ship-launched 
baJlistic  missile  produced  by  an  anomalous  gravity  field.  All  other  error 
sources  are  assumed  to  be  absent.  The  normal  earth  for  which  the  impact 
errors  would  be  zero  is  regarded  as  spherical  and  non-rotating.  The 
standard  atmosphere  is.  accounted,  for.  Both  navigation  and  in-flight  effects 
are  considered,  the  former  under  the  assumption  that  the  ship  remains 
stationary  at  the  launch  point. 


DD  1  JAN ^73  1473  EDITION  OF  I  NOV  BB  IS  OBSOLETE 
S/N  0 102*014*  6601  I 


unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAOE  fSfcMi  Dolt  Mntorod) 


_ UNCLASSIFIED _ 

■  UIHITV  CL  ASSIFICATiON  OF  THIS  PAGEflWnn  Pita  Enftmd) 


The  treatment  is  largely  deterministic,  with  elementary  statistical 
concepts  being  introduced  to  derive  a  description  of  the  impact  errors 
applicable  to  the  case  of  an  unknown  anomalous  field. 

Some  numerical  statistical  resuits  are  presented,  in  the  form  of  a  table 
for  a  particular  trajectory  and  as  graphs  for  trajectories  having  various 
ranges  and  times  of  flight . 


UNCLASSIFIED 


SKCURITY  CLASSIFICATION  OF  TMII 


at/'Vhau  Osls 


